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by resort. to the Differential and Integral Calculus and in practical 
work the Calculus of Finite Differences is almost indispensable. 
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none of them has been written with the special requirements of 
the actuary in view. In beginning his training the student is, 
therefore, confronted by the difficulty of judicious selection and in 
the circumstances it has appeared to the Council of the Institute 
of Actuaries that a mathematical text-book sufficiently compre- 
hensive, with the standard works on Higher Algebra, to provide the 
ground-work of an actuarial education would be of great value. At 
the request of the Council, Mr Alfred Henry has undertaken the 
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AUTHOR'S PREFACE 


ACTUARIAL science is essentially practical in that, whilst it is based 
on the processes of pure mathematics, the object of the worker 
must be to produce a numerical result. 

For this reason it is necessary for considerable prominence to be 
given, in the curriculum of the actuarial student, to the subject of 
Finite Differences, and it thus becomes convenient, in the study of 
those subjects not included under the heading of Algebra, to deal 
with this part of the syllabus first and, notwithstanding certain 
theoretical objections, to treat the fundamental propositions of the 
Differential and the Integral Calculus as being, substantially, 
special cases of similar propositions in Finite Differences. The 
subjects enumerated cover so wide a field that it has been necessary 
to exercise considerable compression and to include only such 
problems as are requisite for a proper knowledge of the subjects 
within the syllabus. 

In the chapter on Probability it will be seen that the numerical 
or “frequency” theory of probability has been adopted. Having 
regard to the practical nature of the actuary’s work, it is thought 
that strict adherence to this aspect of the subject is necessary if 
the student is to acquire sound views from the outset. The subject 
of Inverse Probability has been excluded from the examination 
syllabus in recent years and for this reason it is not introduced into 
the present work. 

In conclusion the author would wish to tender his best thanks 
to many colleagues and other members of the Institute of Actuaries 
for their kind assistance and useful criticisms. In this connection 
he is particularly indebted to Mr G. J. Lidstone, who was good 
enough to read the chapters relating to Finite Differences and 
made many valuable suggestions. 

A. H. 
August 1922. 
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NOTE AS TO SYMBOLS 


In a few cases certain mathematical symbols have been employed, 
the use of which in modern mathematical work may be unfamiliar 
to students. j 


! Used instead of lo to represent a factorial. Thus 
nt=|n=n(n—1)(m—2) «0.1, 


where 7 is positive or integral. 
= Denotes an identity which is true for ail values of the variable. 


! 


(") Used instead of *6,, (n,») or =A 


- Tends to the value ; has the limiting value, e.g. 
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CHAPTER I 


FUNCTIONS. DEFINITION OF CERTAIN TERMS. 
GRAPHICAL REPRESENTATION 


1. When the value of a certain quantity y depends upon, or 
bears a fixed relation to that of another quantity, a, y is said to be 
a function of x, and the relationship is written as y = f(z). 

[Other notations used are uz, vz, p (x), etc.] 

Thus, we may have y=2*, y=a*, y=sing, y=loga. The 
quantities forming the right-hand sides of the equations are all 
functions of «. 

When expressed in this way the relationship of y to « is said to 
be explicit. But if, for example, ax? + 2bay + cy*=0, it is clear that, 
whilst the value of y depends upon that of z, it cannot be determined 
in any case by direct substitution of the value of x. The relation- 
ship in such circumstances is said to be amplicit. In this particular 
example the relationship can, of course, be made explicit by solving 
the equation for y in terms of a. 

The quantities # and y are called variables. 

The quantity 2 is called the independent variable since it may 
assume any value, whereas y is called the dependent variable since 
its value depends upon that of a. 

The independent variable is sometimes referred to as the argu- 
ment. 


2. Functions may, of course, involve more than one variable. 


2 f 
For example we may have y= 42° or w= = ae J. and all these 


variables may be involved implicitly. A function of more than one 
_ variable would be denoted by /'(2, y, 2), Uz.y, ete. 
If we have a function u, such that 
w= ay zt + at y? 2 +... 
and atb+c=a4++c=...=N, 
then w is said to be a homogeneous function of n dimensions. 
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: "Bechibitieu: Co-ordinates. 

3. A function of two variables may be represented graphically 
according to the scheme of the figure shown. 

Let two straight lines of indefinite length OX, OY be drawn at 
right angles to each other. The point O is 
called the origin, OX is called the axis of « 
and OY the azis of y. 

Then if OW is measured along the axis of 
z, equal in value to 2, and at that point a 
perpendicular line PN is drawn equal in 
value to y = f(x), ON is called the abscissa 
and PN the ordinate of the point P. The 
convention is taken that measurements of x in the direction OX are 
considered to be positive and those in the contrary direction nega- 
tive. Similarly measurements of y in the direction OY are treated 
as positive and those in the contrary direction as negative. The 
point P, written for convenience as (a, y), is thus completely deter- 
mined from given values of x and y. 

If for every value of « the corresponding value of y were plotted 
on a diagram such as the above a continuous curve would be 
obtained which would be the graphical representation of the equa- 
tion y= f(x). It is, of course, impossible in practice to plot every 
value of the function, but generally a few values can be filled in so 
as to enable the curve to be drawn by sight. 

It should be noted that the lines OX and OY, and consequently 
the origin O, can be chosen quite arbitrarily and that the position 
of the point P can be fixed, in the manner indicated, with reference 
to any suitable axes of co-ordinates. 


4. The point P can be fixed with reference to its distance from 
two straight lines not at right angles to each other, the distance 
PN = y being measured along a line parallel to the axis of y. 

The values « and y corresponding to a given point P are called 
the rectilinear co-ordinates (or simply the co-ordinates) of the point 
P. Where the axes are at right angles to each other the system 
can be distinguished, if necessary, by referring to # and y as the 
rectangular co-ordinates of P. 

As a rule rectangular co-ordinates are the more convenient to use 
in practice and lead to simpler results. Unless otherwise expressed 
it is to be understood that rectangular co-ordinates are implied. 
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5. The following examples give simple cases of the graphical 
representation of explicit functions. , 


(i) The equation «=a clearly represents a straight line parallel 
to the axis of y and at distance a from it; for the value of # at any 
point is constant and equal to a. 


(ii) The equation y=mzs represents a straight line passing 
through the origin and making an angle 6 with the axis of x, where 


tan =m; since at any point the ratio y: is constant and equal 
to tan @. 


(iii) 


Let AB be any straight line cutting the axes of # and y respec- 
tively at the points A and B, so that OA =a and OB=0. 

Let P be any point on the line AB, of which the co-ordinates are 
(x, y). Then, if perpendiculars PN and PM be dropped upon the 
axes of and y, MP=a2 and NP=y. 


Al a _ON _PB 
<a a OA AB 
y NP _AP 
ang b OB AB’ 
Zy PB+AP 
ob As 
=1. 


Hence “+ =1 is the equation of the straight line AB, 


\1—2 
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6. An implicit function can be similarly represented. For ex- 
ample, it is obvious from the ordinary properties of the circle that 
the implicit relationship «+ y?=a? represents a circle of radius a 
with its centre at the origin. 


Note. The function y=a+ ba+ca*+d2°+... is sometimes 
called a parabolic function, since the equation y=a + ba + ca 
is represented graphically by a curve which is known as a 
parabola. 


7. It does not follow that for every value of « there will always 
be a real value of y. 

Thus, consider the function y? = (a —a)(a—b)(x—c), where 
c>b>a. If x is negative, the right-hand side of the equation is 
negative and y can have no real value. If # is positive and <a, the 
position is the same. If, however, >a and < b, then the right- 
hand side is positive and y has a real value; but when #>b and 
<c, y is again unreal and remains so until « >c when a real value 
of y results for each value of x. 

The form of the curve is shown below, where OA =a, OB=b 
and OC =c. 

In circumstances such as these, 
where one or more parts of a 
curve are isolated from the 
others, the function and the 
curve representing it are saidto 0] 4 BC 
be discontinuous. 


8. It is convenient here: to 
introduce the conception of the 
limiting value of a function, or simply a limit. 

If y =f («) and y continuously tends towards a certain value and 
can be made to differ by as little as we please from that value, by 
assigning a suitable value to z, say a, then f(a) is said to be the 
limiting value of f(a) when a tends to the value a. 

A convenient notation is as follows: 


y—~ f(a) when «a. 


Also f(a) would be expressed as Lt f(z). 
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Thus let y= —- By writing y in the form 1 -" we see that 
by making indefinitely great, we can make the value of y differ 
from unity by as little as we please. 


Thus ee I 


Ct Br oo? 2] YG 


9. We will now give an example of another form of discontinuity, 


# 2 
) shown 
—@ 


and for this purpose we will take the curve y= E 


below. 

Here it will be seen that as «> a, the value of y+ oo. Similarly 
ifx>o, y>l1. 

Thus if we draw two lines, one 
PN parallel to the axis of y and at 


distance a from it, and the other i, Q 
QM parallel to the axis of x and at ¥ 
unit distance from it, the curve will 0 N 


continuously approach these lines 
but will not actually touch them 
except at an infinite distance from the origin. 

Such lines are called asymptotes to the curve. 

In general, actuarial functions are finite and continuous; but in 
mathematical work, as will be seen later, attention to these points 
is necessary in the consideration of certain problems. 


10. Explicit functions involving two variables, and implicit 
functions of three variables, can be expressed in a diagram of three 
dimensions. Thus if we have z=f(a,y) or f(a,y,z)=0 we may 
measure x, y and z by reference to their perpendicular distance, not 
from two lines, but from three planes at right angles to each 
other. 

A useful simile is that of the floor and two adjacent walls of 
a rectangular room of indefinite size. Any point in space is fixed 
by reference to its perpendicular distance from the two walls and 
the floor, 2 and y being the respective distances from the walls, 
and z that from the floor. The conventions as to positive and 
negative values are similar to those previously explained. 
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11. As an example, if z = 2*y?, the values of the function for unit 
intervals in the values of 2 and y are shown in the following table: 


Value of x 


We will conceive the floor as being covered with a linoleum of 
chess-board pattern, the sides of the squares being at unit distance 
apart. The perpendicular distance from the floor of all points in 
these squares is nil and, therefore, z has the value 0. The corners 
of these squares will thus represent the various points (0, 0, 0), 
(0, 1, 0), (1, 1, 0), (1, 0, 0), (2, 1, 0), etc. If we were to erect at each 
corner a peg of height equal to the appropriate figure taken from 
the above table, the tops of the pegs would give points on the surface 
representing the equation z=a*y*, If the value of z were plotted 
for every possible combination of values of 2 and y, we should, 
of course, obtain the continuous surface corresponding to z= ay. 


Polar Co-ordinates. 


12. An alternative method of defining a point in a plane surface 
is as follows. 


Let an origin O be taken and a fixed line OX be drawn from it; 
then the position of any point P is known if the distance OP and 


the angle XOP are given. 
ya J 
‘ 0 


Q 
Thus if OP=rand 2 XOP=6,r and @ are known as the polar 
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co-ordinates of P, and the point P can be written as (r, 0). The 
distance OP is called the radius vector and the angle X OP is called 
the vectorial angle. 

The convention adopted is that the angle XOP is reckoned 
positive if measured from OX in a direction contrary to that in 
which the hands of a clock revolve, and negative if measured in the 
reverse direction. 

Further, the radius vector is considered positive if measured 
from O along a line bounding the vectorial angle, and negative if 
measured in the opposite direction. To illustrate this system, let 
PO be produced to a point Q'such that OQ=OP=r. Then the 
point Q may be written alternatively as (7, 7 + 0) or (—7, @). 


13. The relation between rectangular and polar co-ordinates can 
be easily established. For if OX be taken as the axis of x then 
OY, the axis of y, is perpendicular to it. Also let PN be drawn 
from the point P perpendicular to the axis of a. 


Y 


—— 5 


0 <———- x 


N Xx 
Then, clearly, if x, y be the rectangular co-ordinates of P, 
z= ON = OP cos 6=r cos 8, 
y= PN =O0P sin @=rsin 0. 
Any equation in rectangular co-ordinates can therefore be trans- 


formed into an equation in polar co-ordinates by the above substi- 
tutions. 


14. Three simple examples of the graphical representation of an 
equation in polar co-ordinates are now given. 
(i) The polar equation r=a clearly represents a circle of radius 


a with its centre at the origin; since the radius vector is constant 
_and equal to a. 
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(11) In the diagram shown in § 13, if ON =a and PN be produced 
indefinitely in either direction, then the polar equation of the 
straight line so obtained will be 


rcosé=a, 
since if P be any point in the line 
OP cos 0=r cos 0 =ON =a. 


(iii) Let OA be a diameter of a circle OPA of radius a. Then 
if P be any point on the circle such that OP=r and 2 AOP=6, 
OP =0OA cos 6 = 2a cos 6. 


P 


The polar equation of the circle, if O be the origin, is therefore 
r = 2a cos @. 


CHAPTER II 
FINITE DIFFERENCES. DEFINITIONS 


1. The subject or calculus of Finite Differences deals with the 
changes in the values of a function (the dependent variable) arising 
from finite changes in the value of the independent variable (see 
Chapter I, § 1). 

Many questions arise which can be dealt with on systematic lines, 
but probably the most important problems which require to be 
solved in actual practice, and with which we are concerned at this 
stage of the subject, are the summation of series, and the insertion 
of missing terms in a series of which only certain terms are given. 

It will be convenient to proceed in the first place to some ele- 
mentary conceptions and definitions. 


2. If we have a series consisting of a number of values of a 
function, corresponding to equidistant values of the mdependent 
variable, and from each term of the series we subtract the algebraic 
value of the immediately preceding term, we shall obtain a further 
series of equidistant terms. The process is known as differencing 
the terms of the series, and the terms of the new series are known 
as the first differences of the original terms. By repeating the pro- 
cess with the terms forming the first differences, we shall obtain 
a further series forming the second differences of the original function, 
and so on. Thus if we have f(a) for the first term of the series 
and f(«+h) for the second term, the first difference of /(«) 1s 
S(a+h)—f (a) and is designated Af(x). The second difference of 
F(a) is Af (@ +h) — Af (a) and is designated A?f (x). This may be 
set out as in the following scheme: 


Function First Differences Second Differences - Third Differences 


F - h h) —3f (w+ 2h 
Fw) | ferh)-f(e) — | f(w+2h) yey pees oe te 


a , 4h)-3 3h 
PDN esa a TO) emer es) |e tere, 


f (w@+2h) \’f (e+3h)—f (2h) | f(e+4h)- ve , 
(a+ 
f (w@+3h) | f(wt4h)—f(2+3h) 


10 FINITE DIFFERENCES 


The first term of the series is known as the leading term and the 
terms in the top line of differences are known as the leading dzf- 
ferences of the series. 

It must be clearly understood at the outset that A is merely a 
symbol representing the operation of differencing f(a) once; it is 
in no sense a coefficient by which f(z) is multiplied. This point is 
dealt with again in § 5. 


3. An examination of the character of the series which ultimately 
results from the process of differencing repeatedly, leads to the 
development of certain important theorems. Before proceeding 
further, it will be helpful to give a practical example. 

Example 1. Obtain the differences of the series given by f(a) = 2%, 
where x has all integral values from 1 to 6. 


First Second Third Fourth 
Differences Differences Differences Differences 


It will be observed in the above example that the fourth and, 
therefore, all higher differences are zero; it will be seen later that 
this would equally have been the case had more terms of the series 
been taken. We can therefore construct all the remaining terms 
of the series by a process of continuous addition. 


4. Although most functions with which the actuary has to deal 
are not of the simple character of that shown above, yet it will 
usually be found that the differences of the function for which 
further values are required tend to the value zero and are sus- 
ceptible to treatment by methods which will be developed subse- 
quently. 

The student should obtain confirmation of this fact and insight 
into the character of certain series by taking out the differences of 
tabulated functions such as logarithms, annuity-values, etc. 
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5. Before proceeding to the consideration of the various problems 
which arise, it is necessary to develop certain fundamental formulas. 

In §2 A has already been defined as the symbol of the operation 
by means of which the value of f(x+h)—/(a) is obtained. 
Similarly, it is customary to use the symbol H as representing 
the operation by which the value of f(x) is changed to the value 
f(@ +h), so that 

Ef(a)=f(ethy=f (2) +Af(@). 

It must be carefully remembered that these symbols represent 
operations only and must be interpreted accordingly. Thus £?2* is 
clearly not the equivalent of (#x)?; the former expresses the result 
of operating twice upon the function 2 in the manner indicated 
above, giving a value (a + 2h)?, whereas in the latter case the opera- 
tion is applied once to the function x and the resulting term («+ h) 
is squared. 


6. If, then, these symbols are found to obey the ordinary alge- 
braical laws, they can be dealt with algebraically provided always 
that the results are interpreted symbolically in relation to the 
function which is the subject of the operation. This principle is 
known as that of Separation of Symbols or Calculus of Operations, 

The algebraic laws referred to above comprise: 

(1) The Law of Distribution. 
(2) The Law of Indices, 
(3) The Law of Commutation. 

Taking these laws in succession : 


(1) The symbol A is distributive in its operation, for 


ALA@)+fA@+f@)+-.-J=[LA@thtfh@eth+ferth)t+...] 
—[A@)+h@+f@)+---J 
=[A(et+h)-fi(@)]+[A@t+h)-A@] 
+[fe(at+h)—fs(x)] +... 
=Af,(2)+Afi(~) +Afs(a)+.... 
Similarly the symbol £ is distributive, for 
E (f(a) +fla)+fa(2)+-.-J=[A@+h)tf@th)t+fs@rh)t...] 
=Ef(e)+ Hf, (a)+ Bf, (a) +.... 
(2) The symbol A obeys the law of indices, for in the case of 
positive integers the symbol A™/(«) represents the operation, 
repeated m times, of differencing f(a). 
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Thus A™ f(x)=(AQA...m times) f(z), 
. ATA™ f(z) =(AAA...n times) (AAA... m times) f(z) 
=(AAA...(m+n) times) f(a) 
= Artm f (2), 
Similarly it may be shown that the symbol E obeys the law of 
indices. 
(3) The symbol A is commutative in its operation as regards 
constants, for, if c be a constant, 
A [of(a)]=of(e+h)—cf (2) 
=c[f(e+h)—f(2)] 
=cA f(z). 
The like result can be deduced as regards £. 


7. It follows that, since 
Ef(e)=(1+ A) f(@), 
therefore H=1+A 
and A=£E£-1. 
The two operators are thus connected by a simple relation, which 
will be found later to lead to important results, 


8. As an example of the manner in which the relationship 
between the operations represented by E and A can be utilised in 
the solution of problems, we may take the following: 


Example 2. Prove that 
$0) + af 1+ F fQ)+TS@)+-. 


=¢ | FO) + 2Af +5, MFO) + i, 
Since 
fQ)= #£fO)=(1 + A) f (0); f (2) = £7 f(0) =(1 + AY f(O), ete., 
we have f(0)+ af (1) +5, f@) +5, f(8) + 
= E +a(1+d)+5 (+A) + | 


= [array RY) 
= & [e*] f(0) 


= ¢ E fab tats +...| Ff (0) 


= ¢ | £00) +c Af (0)'+ i A?f (0) + 4 


CHAPTER III 


FINITE DIFFERENCES. GENERAL FORMULAS 
AND SPECIAL CASES 


1. Starting from the relationship proved in Chapter II, §7, it is 
now possible to develop two formulas of the utmost importance. 


2. To express f(a+mh) in terms of f(a) and tts leading dif- 
Serences. 


By definition f (@ + mh) = E™ f(x) 


=(1+ Ayn f(@), 
since by Chapter II, § 7, B=1+A. 


Expanding the above expression by the Binomial Theorem we 
have 


f(et+mh)=[14 mA + (3) Att (F) ass... + (™ abo 
=f ()+ md f (0) + (2) Atf(o)+ (3) Atf(@) +.. 
+ (T) A™S(@) «+. 


Bearing in mind that the symbol A obeys the ordinary algebraic 
laws and that the Binomial Theorem holds for all values of the 
index, positive, fractional or negative, it will be realised that the 
above proof is perfectly general. 

It is instructive, however, to show how the formula may be 
deduced. for fractional and for negative values of m. 


3. Fractional value. Let ~ be a positive proper fraction, and 
let the interval between the given terms of the series be h. It is 


—— to find the value of f (2 + = h) ‘ 
Also let f(at+h)—f(z)=Af (2), 
wd f (e+ =) Faeyr e) 
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In the second case, the unit of differencing has been altered to 
z and, bearing in mind that n is a positive integer, we may write 
at once from the theorem in the preceding article 
f(e+hy=(1+ or f@), 
therefore (1 + A) f(z) =(1 + 8)" fF (2), 
and (1 + AY* f (a) = (1 + 8) f(a). 


Since m is also an integer, it follows that 
(14 5)" f(a) =f (2+ Zh) =(1+d)*F@) 
m 
3 E +e A (=) A? + | F@ 
® 2 
m = 
=f (a)+™ Af(a)+ (=) axv(o)+.... 
4. Negative value. It is desired to find the value of f (a — mh). 


Now, from the preceding theorems, it is clear that 


: (1+ 4)" f(a —mh) =f (2), 
therefore 


f(@—mh)=(1+A-"F@) 
=[1+Cma+("9) a+... J F@ 
= f(a) +(—m)Sf(0)+(~9) AF @ +. 


The above proofs show that the theorem holds universally and 


illustrate how the principle of Separation of Symbols can be applied 
whenever the symbols of operation obey the ordinary laws of algebra. 


5. To express A" f(x) in terms of f(a) and its successive values. 
A™ f(a) = (EH —-1)"f (2) 
=| E"— mE + (3) Be... +(-1)"] f(a) 


=f (w+ mh) —mf(2+m—1h)+ (3) f(@-+m—2h)— vee 


+(—1)™f(a)......(2). 

Alternatively both the above formulas can be easily proved by 
the ordinary methods of induction. 

Formula (1) also follows directly from the ordinary formula of 
Divided Differences (see Chapter VIII). This method has the 
advantage of showing directly the application of formula (1) to cases 
where m has a fractional or a negative value. 
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6. The above formulas are expressed in a form which applies in 
the most general way, i.e. when the interval of differencing is h and 
the leading term is f(#). It 1s clear, however, that by altering the 
unit of measurement the formula will be simplified although the 
result is not affected. Similarly by changing the leading term 
(which process corresponds to shifting the “origin”) so that the 
leading term is expressed as f(0) a further simplification in form 
is made. 

If, therefore, the interval of differencing becomes unity and the 
leading term can be represented by /(0), the first formula can be 
written 

F(m)=Ff (0) + MAFO)+(G) ASO) Foe coves (3). 
An example will make this clear. 

Having given the values of f(10), £(15), £(20), etc. 1t is desired 
to express f(17) in terms of (10) and its leading differences. 

The original formula (1) gives the value of /(10 + 14h), where 
h=5, and therefore we write 


f (10 + 14h) = f (10) +144 (10) es (10) teu, 


where A, A?, ... are taken over the interval h. 

But the same result is secured if the unit of measurement is 
changed from 1 to 5 and if at the same time /(10) is made the 
initial term of the series, for then f(10), (15), £(20), ... can be 
written as F'(0), F(1), F'(2), ... and the required value, viz. f(17), 


becomes F'(1:4) which by formula (3) is equal to 
F(0)+ 140 F (0) +29) A’ F(0)+.... 


7. The above formulas are of general application if sufficient 
terms of the series are known, but it is convenient at this stage to 
consider the particular forms taken by the differences of certain 
special functions. 


8. f(x)=az2", 
Af (a) =a(a+1)"—ax"=a [ na + 4) grt + 1 |. 
The result of differencing has been, therefore, to change the term 


involving the highest power of # from ax” to ana"— (thus reducing 
its degree in 2). 
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Similarly a further process of differencing will reduce the degree 
of « to n—2 and the coefficient of the highest power of « will be 
an(n—1). By repeating the process we arrive at the result that 
the nth difference of az” is independent of x and is equal to a.n!. 
The (n+1)th difference is therefore zero. 


Corollary. It follows that the nth difference of 


ag® + ba + ca® 4+... +k 
is constant and equal to a.n!. 


9. f(a) =2(a@—-1)(e@—-2)...(@—-m+]1). 
This expression is usually denoted by a™, 
A f(#) =(¢4+1) #(a—1)...(a— m+ 2) —#(a—1) (a-2)...(a—m+1) 
= ma (a—1)...(e—m+2) 


= maim), 
Similarly A?f (x) = m (m — 1) a"), 
By repeating the process we arrive at the result 


A™ f(a) =m!, 
which is otherwise obvious from the preceding article since f(z) is 
of the mth degree in 2. 
1 
TO) se+ G2). @ Fm =I)" 
Corresponding to the notation already used, this can be denoted 
by a), 


10. 


il il 


AF@ = Gx @+2)..@4m) 241) @+2)..@tm—l) 
—m 
~ @(@+1)(@+2)...(e+m) 
=— ma\—™), 
Similarly A? ‘i (x) =m (m+ 1) gi-m+2) 
and so on. 
11. Se (z) i a”, 
Af (2) = a" — a? 
=a*(a—1). 


Whence A*f (#) = a* (a — 1)". 


GENERAL FORMULAS AND SPECIAL CASES 17 


12. For many purposes it is convenient to have a table of the 
leading differences of the powers of the natural numbers. These 
can be represented as the differences of [a"],-) and are sometimes 
known as the “ Differences of 0.” 

The following table gives a number of values of the first term 
and leading differences of [A”z"],.5, which, for convenience, can 
be written as A™0”; 


13. A working formula for constructing a table such as the 
above by a continuous process may be obtained as follows: 


An f (0) =f (w +n) —nf (@+n—-1)+(3) f@@tn—2)—.... 


[See wate (2).] 
Therefore, when f(x) = a™, 


a me ee el)” s ( ) (c+n—2)™... 
whence, putting z= 0, 


A"0™ = n™ —n(n—1)™+ (3) (n— 2)" — 
=n[ a" —(n—1)(n— 1+ ("3") (a 2)m—...] 
=n|(l+n—1)"9-(n=1)(140=2)™" 
ag —, (14+n—3)"™> — ~ 


a er aed (4). 
But F(A) =f 0) + AF), 
and A™ f(1) =A" f (0) + A” f (0). 
Therefore 


[Ar a™—2], 2, = [A 2™],2) + [AZ ],25 
= A™10™—1L ep APO"), 
A. TBI, . 
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Hence BOO ae ap Om I). ies (5). 


It follows that the differences of [z™],~) can be constructed from 
those of [z”—],.,, and so on. 
To take an example from the table given above, 
A405 = 4 [A204 + A‘0*] 
= 4,[36 + 24] = 240. 


14. By using the result given in § 9, it is possible to expand 
F(a) in terms of 2, a, af), .... 


Let f(a) = Ap t+ 4,0 + A, + Asa® +... 
Then, putting # = 0, wé see that 
Sr (0) = A,. 


Differencing both sides of the equation, we get 
Af (a) =A, +2A,0+3A,09+4+.... 
Again putting «=0, we find 
Af(0)= A). 
By repeating the above processes, we obtain successively 
A?f(0)=21A,, A®f(0)=31A4,;,... A*f(O)=n!A,, 
whence : 
4,=f0), A=Af0), 4-92 4, ee 

and : 

f(a) =f 0) + af (0) +F Af (0) +... + Af(0) +... (6). 


CHAPTER IV 
FINITE DIFFERENCES. INTERPOLATION 


1. The subject of interpolation is one of the most important in 
Finite Differences and may be enunciated as follows. 

It frequently happens that we have given a number of values of 
jf (#) corresponding to different values of «, and we wish to find 
a value of the function for some other value of a If the form 
of the function is known or can be deduced from the given values, 
the problem is, of course, simple, although in many cases it is more 
convenient to proceed by the methods of Finite Differences. But it 
is frequently the case, especially in actuarial work, that the function 
cannot be expressed, algebraically or otherwise, in any simple form, 
and resort must be had to other devices, 


2. Looked at from the point of view of a problem in graphs, we 
may regard thegiven values of the function as representing a number 
of isolated points on a curve, and it is desired to plot a further point 
corresponding to a given value of the abscissa. 

It follows that if the form of the function (ie. the equation of 
the curve) is unknown, some assumption must be made as to the 
relationship between the different values. The formulas of finite 
differences assume that this relationshipcan be expressed in the form 

y=at ba + cu? + da +... + kam 
This assumes (see Chapter III, § 8) that all orders of differences 
higher than the (n — 1)th vanish, but, as pointed out in’ Chapter II, 
§ 4, this assumption can be made without introducing important 
errors in practically all cases where actuarial functions are involved. 


3. The above equation contains n constants, and therefore n 
values of the function must be known if the values of the constants 
are to be determined. Conversely, if m values only are known and 
the methods of finite differences are to be applied, it must be 
assumed implicitly that all orders of differences higher than the 
(7 — 1)th vanish. | 

4. The most obvious method of procedure is to obtain the n 
equations given by the n values of the function and to find the 
values of the constants therefrom. The assumed form of the function 

2—2 
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is then completely determined and the value corresponding to any 
value of « can be obtained. 

In the majority of cases, however, this is not the most simple 
method of working, for other devices can be adopted which will 
materially shorten the arithmetical work. It is important to note, 
however, that alternative formulas, in which the same values of the 
function are used, lead to identical results. 

In some cases there is scope for the exercise of the ingenuity of 
the solver, but usually the problems fall into the main categories 
which are illustrated in the following examples. 


5. Haample 1. Whenn equidistant values ofa function are given 
and it is required to find the value of some intermediate term or terms. 
This can be done readily by the application of formula (1) of 
Chapter ITI, or by the simpler formula (3). From the given values 
the successive orders of differences are calculated, and the result is 
obtained by direct substitution. 
Thus, taking the numbers living by the H™ table at ages 45, 
50, 55, 60 and 65, it is required to find the value for age 57. 
In conformance with formula (3) the given values can be denoted 
by f(0), f(1), ..., so that the required value is /(2°4). Then 
2: Ax x 1-4 


Ff (24) =f(0) + 24 AF (0) + A’ f(0) 
+ There® 4 (0) 4 ATER KF sory 


The working is as follows: 


A? f (2) A* f (2) Af f (z) 


— 1106 — 389 +75 
— 1495 —314 
~ 1809 


From above: 
f (24) = 77918 + 2°4 (— 5123) + 1°68 (— 1106) 
+ ‘224 (— 389) — 0336 (+ 75) 
= 77918 — 12295°2 — 18581 — 87:1 —2'5 
= 63675'1. 
The value given by the table is 63677. 
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The difference between the interpolated value and the true value 
is due to the fact that the interpolation curve, which is based on 
the assumption that all differences of higher order than the fourth 
vanish, represents only approximately the true function. 


6. Haample 2. When the values given and the value sought 
constitute a series of equidistant terms. 

If there are m terms given of which n —1 are known, then, as 
explained in § 3, it must be assumed that the (n —1)th order of 
differences is zero. 

Thus, using formula (2) of Chapter IIT, we have 


Ant f(0)=0=f(n—1)- (n—1) f(m—2)+ ("9") F(n—-8)--.. 
+(-1"f). 


In this equation there is only one unknown quantity and its 
value can, therefore, be readily obtained. 
For example, if 
F (0) = log 3°50 = 54407, 
Ff (1) = log 3°51 = *54531, 
J (2)= log 3°52 = 54654, 
FJ (4) = log 3°54 = 54900, 
and it is required to find log 3°53, 1.e. £ (3). 
From above: 
A‘ f(0)=0=f (4) — 4f (8) +6 (2) - 41) +40), 
F (8) _f (A) aL i 4f (1) +f (0) 


= 54777, 
which agrees with the true value to five decimal places. 


whence 


7. Eaample 3. If more than one term is missing from the com- 
plete series, a somewhat similar process may be followed. Thus, if 
two terms are missing, only (n—2) terms are known and the (n—2)th 
order of differences must be assumed to vanish. It is then possible 
to construct two equations: 

Ar f (0) =f (n—2)—(n—2) f(n—3) +... +(—1)"? f (0) = 9, 
A®? f(1) = f(n—1)—(n— 2) f(n —2) +... +(— 1)? fC) =0. 
From these equations, the values of the two unknowns can be 

calculated. 

Similarly if a larger number of terms is missing, the method 
can be extended. 
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8. Haample 4. If several equidistant values are given, together 
with one isolated term. 

For instance, if three values f(0), f(1) and f(2) are given, to- 
gether with a further value f(/). Having four values of the function 
it must be assumed that the fourth order of differences is zero and 
it remains to find the values of the other three leading differences. 
The first two leading differences are obtained at once by differencing 
the first three terms of the series, and the value of the third dif- 
ference is then given by the equation 


f(t) =f (0) + hAF (0) + (3) A*F(0) + (5) A*F(O). 
For example, taking the numbers living by the H™ table at ages 
45, 46, 47 and 50, it is required to find values for ages 48 and 49. 
A A? 
f(0)=77918 —-954 —32 
f(1)=76964 —986 
f (2) = 75978 
f (5) = 72795 
f (5) =f (0) + 54 f (0) + 10A?F (0) + 10A*F (0), 


arf (0) = £EV=LFO)+ 54 f 0) + 108° 0) 


_ 72795 —[77918 — 4770 — 320] 
= 10 | 
=~ 33, 


The table is then completed by addition. Thus: 


F (2) 


77918 
76964 
75978 

74956°7 
73896°8 
72795°0 


The work is checked by the reproduction of the value for age 50. 
The tabular values for ages 48 and 49 are 74957 and 73896 


respectively. 
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Values precisely the same as those obtained above would have 
been given if the two missing terms had been inserted by the 
method described in Example 3. It is instructive to confirm this 
by actual calculation and to compare the two methods of procedure. 


9. Example 5. Subdivision of Intervals. 


This problem arises when a series of equidistant terms of a series 
is given (usually every fifth term or every tenth term) and it is 
desired to find by interpolation the values of all the intermediate 
terms. 

The simplest method of procedure is to calculate from the given 
values the differences corresponding to the individual terms of the 
series (the subdivided differences) and thence to construct the table 
by summation. The calculation is checked by the reproduction of 
the values of the original terms. 


Thus assume that the given terms are f(0), f(1),... f(5) and it 


is desired to complete the series f(0), f(4), f(2), etc. It is con- 
venient to adopt the notation 


f()-f (0) =F (0), 
and f(t) —f (0) =8f (0). 
The problem then becomes to express of (0), &f (0), ... in terms 
of Af (0), A?f(0), .... 
Writing f (1), f(2), ... in terms of the subdivided differences, 


0)= 
Ae Fo} 58f(0)+ 1082f(0)+ 1083f(0)+ 58*f(0)+ 35 f (0) 
f (2)="F(0) +108 (0) + 458%f(0)+ 1208°F(0)+ 2108*f(0)+ 25288F (0) 
f (3) = F (0) + 158f (0) +10582F(0)4+ 4558°F(0)+ 136584f(0)+ 30030°F (0 


f (4) =f (0) + 208f (0) + 19082F (0) + 114033f (0) + 484534F (0) + 1550485 to} 

f (5) ="f (0) +258 (0) +30032f (0) + 23003%F (0) + 1265084F (0) + 531308°F (0) 
Differencing successively both sides of the equation, we have 

(0)= 58f (0 1082 (0 108°f(0)+  584f(0)+ 65f (0 

r= parie asrt0)4 nea eae at aa es13ey 0) 
Af(2)=  B5Sf(0)+ 608°F(0)+ 3353°7(0)+115534f (0)+ 27518°F (0) 
Af(3)=  5df(0)+ 8582f(0)-+ 6858%/ (0) +348084f (0) + 125018°F (0) 
af(4)=  5df(0)+ 1108%f (0) +11608%f (0) + 78053*f (0) + 376269°F (0) 
2f(0)= 2567F(0 1008°f (0)+ 200847 (0)+ 25085/(0 
erie ero: paper (0) « ge Pc 
A*f(2)—= 2582F(0)+ 3508%F (0) + 232584F (0) + 975085 (0) 
A*f(3)= 2588 f(0)+ 47538 (0) + 432594f (0) +-251258° f (0) 
Af (0)= 12533f(0)+ 75034f(0)+ 225085 (0) 

A®f(1)= 1258°f (0) +13759¢f (0)-+ 72503°F (0) 

A®f(2)= 1258%f (0) + 200084f (0) + 1537585 Ff (0) 

A‘f(0)= 62534f (0) +-500085f (0) 

A‘f(1)= 6258¢f (0) +81250%f (0) 

AS f (0) = 312585 f (0) 
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Whence the values of 6/(0), 5°f(0), ... 6°f(0) can readily be 
obtained. 


10. Alternatively the formulas for 6, 6?,... can easily be written 
down by using the method of Separation of Symbols. 


Thus (1+8)f(#)=(1 +A)/ (2). 
Therefore (1+8) f(#)=(1+A)* f (a), 


and 8 f (a) =[(1+4)* — 1] f (@) 
= [2A —08A? + 048 A8...] f (2). 
Hence & f (x2) = [2A — 08 A? + 048A3...P f(x) 


» = [-04A?—-032A8 + 025644...) f(2) 
and so on. 
For convenience the coefficients of A, A*,... occurring in the 
values of 6, 6*,... are given, for the intervals 5 and 10, in the fol- 
lowing tables. 


Subdivision into 5 intervals 


Coefficient of 
Value of 


+ °025536 
— 02112 
+ °00960 
— ‘00256 
+ ‘00032 


Subdivision into 10 intervals 


Coefficient of 
Value of 


A’ Aé A® 


+0285 | — °0206625 + °01611675 
— ‘009 + °007725 — 0066975 
+°001 — 00135 + ‘0014625 
+ 0001 — 00018 
+ ‘00001 
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11. The following example gives an illustration of the method 
of working. 

Given the present values, at 3 per cent. interest, of an annuity of 
1 per annum for 20, 25, ... 45 years, it is required to find the inter- 
vening values. : 

We have 


F (2) 


14:8775 
17°4131 
19°6004 
21°4872 
23°1148 
24:5187 


Applying the factors given in the former of the above tables, we 
have 
df (20) = [2A — 08A? + -048A% — -0336A‘ + °025536A5] f (20) 
= ‘2 x 2°5356 + 08 x ‘3483 + 048 x ‘0478 + ‘0336 x ‘0065 
+ 025536 x ‘0007 


= 5375146752. 

Similarly 
6? = — ‘015642784, 64 = — 000012192, 
&= 000451520, &= 000000224. 


The table is then constructed by addition from these leading 
differences as shown below. 

It is necessary to consider how many decimal places should be 
retained in the working. This depends, in the first place, on the 
degree of accuracy desired in the result. Thus if four decimal 
places are required in f(a), our result must be given to at least 
one place more. Further, the range of the formula is 25 terms, and 
in the final term the values of 6, &,... 6° are multiplied respectively 
by (7 ) ; s. Ss (@), which are respectively equal to 25, 300, 
2300, 12650, 53130. In view of the magnitude of these coefficients, 
we shall need to retain at least six decimal places more in the value 
of 6° than are required in the final value of f(). Since the co- 
efficients of the lower orders of differences are smaller than 53130, 
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a correspondingly smaller number of decimal places can be retained. 
In practice, however, there is little to be gained by cutting down, 
unless only a rough result is required. 

The working of the first five terms in the example is shown 
below. It will be noticed that the accuracy of the work up to this 
stage 1s checked by the exact reproduction of the value of f(25). 
The interpolated values agree exactly with the true values to four 
places of decimals. 


14°8775 


15°415015 | 5218719 
15°936887 | 5066806 
16°443568 | °4919287 
16-935497 | -4776041 
17-413101 | °4636951 


12. Example 6. Lagrange’s Theorem. 


We have now to consider the construction of an interpolation 
formula which will apply when the given terms of the series are 
not equidistant. 

Let n values of the function be given, namely 


F(@), FO), FO), «FM. 


Then the function must be assumed to be a parabolic function of 
x of degree (n—1). (See § 3.) Assume therefore that the function 
can be represented as 


f(z)= A(a—b)(a@—-c)...(e@-n) 
+ B(a-—a)(«—c)...(«@—n) 
+ C(x-—a)(@—b)...(a@—n) 
+ etc., 
there being n terms in all, each composed of (n —1) factors multi- 
plied by a constant, the values of the constants having yet to be 
determined. It is clear that the right-hand side of the equation 
is of degree (n — 1). 
To find the values of the constants we proceed as follows: 
Put = a. 
Then f(a)=A (a—b)(a—c)...(a—n). 
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f (4) 


Therefore ?GanGee. AG» : 
Sl J (6) 

Similarly B= @ ~ayb—c)..(b— wy’ 

and so on. 


Substituting these values of A, B,... in the peas equation 
we have 


(2 —b)(w- (2 —a)(a—0)...(@—n) 
I (2) = f(a eB e=9-69 +70 4 ar 
(w7—a)(a@—6)(x—¢).. 
nie aE a)(n—b)(n—c).. AP 
By an obvious transformation, the formula can be put in a some- 
what simpler form for calculations, namely 


f (2) = (2) 
(c—a)(x—b)...(a—n) (a—a)(a—b)(a—c)...(a—n) 


f(6) f(r) 
= 6 a)(#@—6)...(b— ny ne a)(n—b)(n—c)...(@—- i 


In memorising the formula it should be noted that the denomi- 
nators are made up of the product of the algebraic differences of 
the values of the variable, the term (a — a) being replaced by (x — a) 
and so on. 


13. The formula is somewhat laborious to apply, and careful 
attention to signs is required, but it is convenient to use where only 
one or two unknown values of the function are required. Since the 
assumptions underlying it are precisely similar to those previously 
explained, its use in any particular case will give identical results 
with those which can be obtained by the use of the ordinary 
methods of Finite Differences where a sufficiently high order of 
differences has been taken into account. 

To illustrate this point and to provide an example of the use 
of the formula, we will calculate by Lagrange’s formula the value 
for age 49 in Example 4. In this case we have 


f(0)=77918, f(1)=76964, f(2)= 75978, f(5)=72795, 
and it is required to find the value of f(4). 
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We have accordingly 
f (4) f) 
G—0) @= 14-3) 4-5) G0) 0-1) (0-2) 0-5) 


Se) Or a ee 
(1—0)(4-—1)(1—2) (1-5) | (2-0) (2-1) (4— 2) (2—5) 
. f(6) 
(5 —0)(5—1)(5 —2)(4—5)’ ; 
or, a S(4)=- ASO) +f — SO) — dof (5) 
Whence we find f (4) = 73896°8, as before. 

It should be noted, as a check on the formula, that the sum of 
the coefficients of the terms on the right-hand side of the equation 
must equal the coefficient of the term on the left-hand side of the 
equation. 


14. Problems of interpolation between terms at unequal inter- 
vals can also be dealt with in a simple way by the formulas of 
Divided Differences (see Chapter VIII). 


CHAPTER V 


FINITE DIFFERENCES. CENTRAL DIFFERENCES 


1. It has already been stated that in interpolating between given 
values of a function the form of the expression connecting these 
values is assumed to be parabolic, and that this assumption is usually 
only an approximation to the truth. It remains therefore to be 
considered by what methods the best result can be obtained by the 
processes of Finite Differences. 


2. In developing the formulas of this chapter, it will be assumed 
that a number of equidistant values of the function are given. 
Let us assume further that it is desired to interpolate a value f(x) 
intermediate between f(0) and f(1). It is clear that our knowledge 
of the shape of the curve on which the points lie is increased if we 
are given values of the function lying on both sides of f(0), and 
that generally the best value of f(z) will be obtained, if a limited 
number of terms is to be used, when the required value occupies 
as nearly as possible a central position in regard to the terms used 
in the interpolation. 

The formulas of Central Differences are designed to give effect 
to these considerations. 


3. The more familiar formulas of Central Differences are as 
follows: 


STIRLING’S : 
f(a) =f(0) +a a 5 Atf(-1) 


Fie yapensayen, Pe) arf(-2) 


oe Lie 4) A® f(— Dalle 3) 


A (1) 
ot ne 
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BESSEL’S : 
f@=LOV® 5 ear) 
a — 1) A?f(- as A? f(0) 4 (x — — —1) A+ f(—1) 
(w+1)2(w-1)(w—2) A‘f(-1) + A*f(—2) 
=< a 
“ a A’ f(—2) 


(a+ 2)(a+1)2(x—-1)(x—2)(a-3) ASf(—2) + AP(-8) | 
J 6! 2 


GAUSS’: 
Sf (@)=f(0) + cAf(0) + A*f(-1) + 
¢ CED RE=Ml— 9) ar¢_9y 


A*f(—1) 


x2(a—1) (7+ 1)a(x—-1) 
2! 3! 


ee 2) 4... (8). 


EVERETT’: 
F (@) = af (1) + —— A*f (0) + 
HEDENE Pa ayy., 


a 1) a (a 


—1)(#-4 
7e aE ayy) 


4 yf (0) + y _—_ 1) A?f(— 1) _" yyy’) Atf(— 2) 
> | 24, 29 
VY es Vy ) asf(—8) + 
[where y= 1—<]. 


4. These formulas can be obtained in various ways from the 
ordinary formulas of advancing differences. Once, however, the 
scheme of differences entering into a formula is settled, the co- 
efficients can readily be calculated by the method of Separation 
of Symbols. An example may be given of the demonstration of 
Gauss’ formula by this method. 
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Eaample 1. 
To express f(z) in terms of f(0), Af (0), A?/(—1), A®f(—1), .... 


Let (2) = A, f (0) + A, Af (0) + A,A?f(—1) +... 


Then, since 
arf (-1)=75- 4 fs AY4(-D= SFO); ete, 


A? DN 
a = a 
Aw 


G+aya*4 


A? 
es "(ay 


Multiplying up by (1+ A)™, and equating coefficients of A”; 


4” _(rt0-Dirta—2).. .(#—r+l1) 
a (2r—1)! : 


And, multiplying up by (1 + A)’, and equating coefficients of A”, 


1 —r+l1 
a (r+2)(r+a—- = ya (@—rt Dy 
Hence, by subtraction, 


_rte-lCte—2)...@—7) 
2r! : 


Therefore f (a) =f (0) + zAf(0) 
— Af (- 1) + SEDEEWY pep) 4... 


The other formulas should be proved, in a similar way, as exercises 
by the student*, 


5. The formulas of central differences, although in a different 
form, are intimately associated with those of advancing differences. 
For example, if an interpolated value is calculated by using the 
first three terms of Stirling’s formula, it is obvious that the values 


* See J.I. A. Vol. 50, pp. 28-33, 
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of f(—1), f(0) and f(1) are brought into the calculation. It is 
easy to show that the result is identical with that obtained by 
using the first three terms of the advancing difference formula 
starting with the term f(— 1). 

It may be observed that the first two terms of Stirling’s formula 
also involve three values of the function; the third term merely 
introduces the correction necessary to make the formula true to 
the order of differences (i.e. the second) implied by the use of three 
terms of the series. Thus, as the (2r)th and the (27 + 1)th terms 
of Stirling’s formula both involve the use of (27+ 1) values of 
the function, there is ordinarily little advantage in using the extra 
(2r + 1)th term in any calculation. 

Similarly in Bessel’s formula no material increase in precision is 
gained by using 2r terms rather than 27 —1 terms. 

Gauss’ and Everett’s formulas are each true to the order of 
differences involved and for general use they would appear to be 
the best of those propounded. 


6. In view of the remarks at the beginning of the foregoing 
article, it may well be asked what are the advantages of central 
difference formulas, as compared with advancing difference formulas 
so chosen as to make the interpolated term as nearly as possible 
the central term of those employed. It may at once be said that 
the theoretical advantages are small but that the practical ad- 
vantages may be considerable. Thus if it be desired to introduce 
further terms of the original series into the calculation, the 
original calculations relating to the central difference formulas 
hold good, and the values of fresh terms of the formula can be 
calculated until the desired degree of approximation is attained. 
If however an advancing difference formula is used, the introduc- 
tion of fresh terms of the original series, while retaining the 
interpolated term in a central position, necessitates the changing 
of the origin and the recalculation of all the terms of the formula. 
An example will make this point clear. 


7. Example 2. Required to interpolate the value of a unit 
accumulated for 17 years with compound interest at 5 per cent. 
per annum, having given the values for 0, 5, 10,... 30 years. 

For central difference formulas we must take our origin at 15 
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years, and we will take 5 years as the unit. Thus we get the 
following scheme : 


: 02110 | 00580 
1:27628 | ° . 02690 | 00745 


1°62889 
2°07893 
2°65330 
3°38635 
4°32194 


As the unit of time is 5 years, the required value is repre- 
sented by f(°4). 


Central Differences. We will use Gauss’ formula, i.e. 
f(2)=f 0) + Af (0) 
a(@#—1 z+l)ae(x-—1 
e ( 1 ) ar¢(—1) +§ ld ) Av f(—1) 4... 
The successive terms, with the corresponding values of f(:4), are 


Ff(0)= 2:07898, Ist approx. = 2:07893 
‘4Af(0)= °‘22975,2nd ,, =2°30868 


2x— atf(-1)= — 01492, 8rd, = 229876 

WAX ERS asf) =—00192,4th , = 229184 

a ae A‘f(-2)= 00017,5th , =229901 
cia SAC TE 2)= -00002,6th , =2-29203 


It is clear that no further terms would affect the calculated 
value. The true value 1s 2°29202. 


H. T.B.1. 3 
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Advancing Differences. 
1st approximation f(0) = 2:07893, 
2nd approximation f(0)+ ‘44 f(0) = 2:30868, 


3rd approximation 
‘4 x *4 
f(-I)+14af(-1) +! > 


4th approximation 


A*f(—1) = 229376, 
Al-I+14af(-1) +S 4 are 


——— A*f(—1) = 229184, 


++ 


5th approximation 
2°4 x 1°4 
f(—2) + 240 f(— 2) += 2 — a f(— 2) 


9:4, ° ° 
a xa 4 As f(—2) 


p SERRE ERAS asp (—2) = 229200, 


6th approximation 
2°4x 14% °4 x—-"6 x —-16 


(5th approximation) + 120 


A®f(—2)=2°29202. 
It will be observed that in proceeding to the 8rd and 5th ap- 
proximations using advancing differences every term in the formula 
has to be recalculated, whereas, in the application of the central 
difference formula, terms already calculated hold good whatever 
be the degree of approximation. 

It should be noted, however, that both formulas give mathe- 
matically the same results, the difference of a unit in the final 
figure being due to the use of only five places of decimals 
throughout. 


8. As regards other practical points, it may be observed that 
the numerical coefficients in the central difference formulas are 
smaller than those in the advancing difference formulas (see 
Example 2). 

Other advantages arise in special cases. Thus Bessel’s formula 
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can conveniently be applied for the bisection of an interval, since 
the alternate terms vanish, giving 


_f(0) +f) : A*t(- 1 AS, 


fHQatOR _- 
aisle 2)+Af(-1) , 
128 2 _ oF 


Everett's formula gives the same value. 


9. It should be noted as regards Everett’s formula, that in cal- 
culating a series of values the work is nearly halved since it will 
be found that terms in the formula can be made to do duty twice, 
“‘¢” terms reappearing in the calculation as “y” terms. 

This will be seen at ae for, 


i ofay+2 » arf(0) + FE DEH argc 1) tess 


+yf(0)+ Terea arf(-1) + LPO Darga)... 
and OS 
flt+y=yf(2) oo area) PD) azo) .. 


x(a? —1) (a? — 
5! 


+ af (1) +74 —— f(0)+ 4) Atf(—1) 4.0, 


the last line being end with the first. Thus, if we are inserting 
terms in a series by subdividing the interval into five equal parts, 
x='2,'4,... and y=‘8, ‘6,.... Therefore half of the terms used 
in the calculation of f(-2) can be made to do duty in the calculation 
of f(1'8), and similarly for the other terms. 


10. An example will indicate the method of working. 


Example 3. Using Everett’s formula, interpolate the missing 
terms in the following series, between f (40) and (50). 


3—2 


8 


_ 
— 


2 
°4 
6 
8 
2 
*4 
6 
8 
2 
°4 
6 
8 
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The coefficients of the several terms in Everett’s formula are 


2 — 032 006336 
“4, — 056 ‘010752 
6 — 064 011648 
8 — °048 008064 


The work may be arranged in tabular form: 


2x? —A)lr?_4 Sum of first 
pe A‘f(—1)] three terms diene 


| A’f (0) r 
2)+(8)+(4)| forme | (5) +(6) 

(5) (6) (7) 
197°6 
395'8 
5953 
7968 


241°0 
483:0 
726°8 
9735 


308°1 
617°7 
930°5 
1248°1 


_ Columns (2), (8) and (4), which represent the first three terms of 
the formula, are obtained by ordinary multiplication. Column (5) 
gives the sum of these terms. From what has been said above, it 
is clear that column (6), which represents the second set of three 
terms of the formula to fourth central differences, is obtained by 
writing down, in reverse order, the values of column (5) applicable 
to the previous group of terms. The addition of columns (5) and 
(6) then gives the desired result. 

The given values of f(#) have been taken from the tabulated 
values of the probability of dying in a given year of age according 
to the H™ mortality table, multiphed by 10° 

The tabular values for the interpolated terms are 1038, 1081, 
1122, 1172, 1281, 1845, 1415, 1490. The small differences between 
these values and the interpolated values are due to the fact that 
the H™ table was constructed by means of a mathematical formula 
which is only approximately represented by Everett's formula. 
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11. Another method of applying the principles of central 
differences is to express the required function in terms of known 
values of the function among which it occupies a central position. 
This can conveniently be done by Lagrange’s formula. The for- 
mulas are of two types according as the number of terms involved 
is odd or even. Thus we have by Lagrange: 


Number of terms 2n +1. 
3-term formula, 


ae (a — Lin 2(a+1) , “ee — 1) oe ogee 


5-term formula, 


fe) fC) FED, fO_ f@ , FO 
x (a’?—1)(a—-4) 24 (+2) 6(@+1) 4c 6(a-1) 24(a@-2) 


7-term formula, 


J (#) f(= 8) Mme) « See 


ee 


w (a — 1) (a#—4)(a—9) 720(@+3) 120(a+2)” 48(a+4 1) 


Jo Oe oe 
36a * 98 @—1) 1) 120(¢—2) ' 720(@33)"°"" j 


Number of terms 2n. 


4-term formula, 


fe) _§_ f-9,f-D)_ f@ , f® 
@-DH@-H) 6 @+H *2@th We-H* 6@-H 


6-term formula, 


f(@) A Ee 


oo —_ 


@-D@-)e@—B ~W0@+H* 4@+H 


Le « WD aa B®) fo 
~T20+h 1 B@-) Be_-p* D@—-H 


12. These formulas, of course, yield identically the same results 
as other central difference formulas embracing the same terms. 
To illustrate this we will recalculate the value of f(‘4) in the 
example given in § 7. 
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Example 4. See Example 2. Seven terms are given, the formula 
will therefore be 


f(A) _f(-3) _ f(-) , fCW 


‘4-16 —1)(16—4)(16—9) 720x34 120x24° 48x14 


FO), XY F (2) F(3) 


~ 36x47 48x—-6 120x-16* 720x—26 


or f (4) = — 0046592 f(— 3) + 0396032 f(— 2) — 169728 f(—1) 
+ ‘792064 f (0) + 396032 f(1) — 0594048 (2) + 0060928 /(3) 
=+ -05055 — 00466 
164665 :27647 
105079 -20117 
‘02633 
+ 2°77432 — °48230 
= 2°29202 as before. 
Note, as a check, that the algebraic sum of the coefficients of the 
terms on the right-hand side of the above equation is unity. 


13. For the sake of completeness it is necessary to refer to a 
system of notation in connection with central differences which 
was introduced by W. S. B. Woolhouse and is still in use to some 
extent. This system of notation is compared with that used in the 
previous chapters in the following scheme: 


Ordinary Notation Woolhouse’s Notation 

f(-2) i-@) 

Af(-2) a2 
om) A’f(—2) i b.1 

Sf (-1) A?f(— 2) ay Cty 
F (0) Af (-1) Atf(—2) F(0) (a) b (G0) % 

Af (0) A*f(—1) ay Cy 
FQ) af (1) A*f (0) F(2) b, 
Ff (2) F (2) 


where A f(— 2), A? f(— 2), ete. are denoted by a_,, b_,, etc. and 
Ay = $4. + 24,1), 
= $(c_,+C4,). 
Under this notation Stirling’s formula to fourth differences is 
xg? —1 x (a? —1 
I (@) =f (0) +20, +5, +29) 31 ) Cot ag, 


=f(0)+ (a - s) a+ (54 — Ott at + oh ...(11). 
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Similarly Gauss’ formula can be written 
Pit= 70) + ait, + = Sy (a+ = 7 


4 et ate De =— (12). 


14. Another system of notation, which is extensively used, is 
that due to W. F. Sheppard. Two operators 6 and » are used, such 
that 


Sf(-H)=fO-fE-Y, 2fH=3(/O+/O)) 
8f(G)=f(1)-FO), wd F(0) = $ [87 G) + Of (— 9), 
etc. etc. 


This notation, although somewhat complicated, gives the usual 
central difference formulas in very convenient forms. 


CHAPTER VI 
FINITE DIFFERENOES. INVERSE INTERPOLATION 


1. In direct interpolation a series of values of the function is 
given and the problem is to find the value of the function corre- 
sponding to some intermediate value of the argument. 

In Inverse Interpolation the problem is reversed and it is required 
to find the value of the argument which corresponds to some value 
of the function, intermediate between two tabulated values. 


2. In certain cases of mathematical functions the desired result 
can be obtained by direct calculation. Thus if 


y =f (2) =a", 
gamer 
log a 


and the value of x can be found equivalent to any given value 
of y. 

Where this is not the case various methods can be adopted. 
These will be examined in order. 


3. Let y= f(x) be the given value. Then 


a (a — 


y=f@) =f) + 2d f0)4 “FD wp) +... 


If it be assumed that the higher orders of differences vanish, and 
that the values of A, A?, etc. are obtained from the given terms of 
the series, then we have an equation in x which can be solved by 
the usual methods. 

The disadvantages of this plan are firstly that an equation of 
higher degree than the second is troublesome to solve, and secondly 
that for certain functions the degree of approximation may not 
be very close. Since a quadratic equation employs only three terms 
of the series, it often happens that no close approximation can be 
obtained. In all cases the intervals between terms should be as 
narrow as possible, so that accuracy may be increased and the use 
of higher orders of differences obviated as far as possible. 


INVERSE INTERPOLATION 41 


4. This difficulty of solving an equation in # of higher degree 
than the second can be overcome in two ways. Assume, for purposes 
of illustration, that four values of the function are given, viz. f(0), 


F (1), f (2) and f (3). 
Then 


F(@) =f (0) + ©AF(0) +7 F—9 ary) 4 °F DE? aro) 


No further differences can be calculated and therefore, since f(x) 
is known, the corresponding value of x is found by the solution of 
a cubic equation in # The solution of the cubic can however be 
avoided by proceeding as follows: 

Taking three terms only at a time 


F(a) =f (0) +20) + =F— arg (0), 
and f(a) =/(1) + @-1) Aft) + P22) rg), 


The third difference error in the first equation is 
w(a—1)(@—2 
(a 7 ) A*f(0), 


and, in the second equation, 
(«—1)(@—2)(e—3) 
sea a 


If now both sides of these equations be multiplied respectively by 
(3 — x’) and a (where «’ is a rough approximation to the required 
value, obtained by inspection) and the equations so weighted be 
added together, a new quadratic equation in « will be formed from 
which the third difference error will be practically eliminated. The 
work of solving a cubic equation has been avoided, but all terms 
have been used without sensible loss in accuracy. 

If the mere arithmetic mean of the equations were taken, with- 
out weighting as above, it is possible that, in certain cases, a worse 
result would be obtained by taking four terms instead of three. 


5. Alternatively the solution of the equation may be obtained 
by successive degrees of approximation. 

Thus, taking the above equation and neglecting differences of 
the second and higher orders, we obtain as a first approximation 


the value 2,, where 
he (1). 
Af(0) © 


i= 
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A further approximation is obtained by taking second differences 
into account and writing 2, in place of x in the equation, thus giving 
er oo. a @) 

AO) Fi FO) 2). 
When third differences are taken into account 2, is written for x, 
giving 

nS age pe TE. tee ee (BY. 

Af (0) +4 (2-1) A*F(0) + 3 @2— 1) (a — 2) A8F(0) 

These processes can be repeated until the desired degree of 
approximation is reached. The method has the disadvantage of 
being somewhat laborious. On the other hand it has the advantage 
that an error of calculation at an early stage does not vitiate the 
result, being rectified by the further approximations. 


Le 


6. A different method of procedure is to treat # as a function of 

J (a). Thus since 
y =f (2), 
we may write c= p(y). 

We therefore treat « as a function of y and, since the given 
values of y (i.e. f(0), f(1), etc.) will usually represent unequal 
intervals of the variable y, we must resort to interpolation by such 
a method as Lagrange’s or Divided Differences, in order to obtain 
our value of « (i.e. ¢ (y)) corresponding to the given value of f(z). 


7. The following example is worked out in each of the above 
ways. 
Example. Find the number of which the log is 2}, having given 
log 200 = 2°30103 


210 = 2°32222 
220 = 2:°34242 
230 = 2:36173. 


Method I. 


2°30103 


2'32222 
2°34242 
2°36173 
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f(x) = 233338 = 230108 + 021192 — 00099 ~ me aa 1) 


a (a — 1) (a — 2) 
+ ‘00010 6 


Or, by reduction, 
x? — 32°7 a? + 1303'1 2 — 1938 = 0. 
Whence, solving the cubic, 
x = 15443. 
Since the initial value of 2 is 200 and the unit of measurement 


is 10, the result of the calculation is to give 215°443 as the required 
value of a. 


8. Method II. 


f(a) = 230108 +-02119%— 00099 =») 


also, f(a) = 2°32222 + ‘02020 (a — 1) — 00089 o@-De=9 . 


The first approximation to the value of w# is # =1'5, so that 
3—a'=1°'5. Since the values of 3—' and 2’ are approximately 
equal, we may take for our “weighted” equation the arithmetic 
mean of the above equations, giving 


fF (&) = 233333 = 2°30108 + 02161 « — 00047 x. 
Whence, solving the quadratic, «= 15443 as before. 


9. Method III. 


[st approximation 
2°33333 —2°30103 |. 
“= — (ao | 5243, 
2nd approximation 
2'33333 — 230103 


202119 — 4 x 52483 x 00099 


3rd approximation 


= 1°5432, 


- 233333 — 2:30103 
** 02119 — } x *5432 x 00099 + 4 x 5432 x (5432 — 1) x 00010 
= 15442, 


which differs only slightly from the value obtained by Methods I 
and IT. 
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10. Method IV. Under this method we may consider the data 
to be as follows: 


2°30103 
2°32222 
2°34242 
2°36173 


It is required to find the value of f(2'33333). By Lagrange, 
(01111) (— 00909) (— 02840) x 200 
(— 02119) (— 04139) (— 06070) 
i (03230) (— 00909) (— :02840) x 210 
(02119) (— °02020) (— :03951) 
(03230) (01111) (— 02840) x 920 
(04139) (02020) (— 01931) 
(03230) (01111) (— 00909) x 230 
(06070) (03951) (01931) 
= — 10°7748 + 103°5416 + 1388764 — 16°2006 
= 215°443, 

This result agrees with the values previously found. The true 
value is 215°442 and it will be seen, therefore, that all of these 
methods give results in this case which are closely approximate to 
the true value. 


f (233338) = 


11. Methods I—III differ from Method IV, in that the first 
three methods assume that f(z) is a parabolic function of z, whereas 
Method IV assumes that # can be expressed as a parabolic function 
of f(z). Both these assumptions may be sufficiently accurate in 
many cases, but in other cases an inspection of the trend of the 
given values of the function may indicate which assumption is to 
be preferred. 


12. By the use of the formulas of Divided Differences given in 
Chapter VIII, Methods I—III can be applied to cases where the 
given values of the function are at unequal intervals. 

The application of Method IV is, of course, perfectly general. 


e 


CHAPTER VII 
FINITE DIFFERENCES. SUMMATION OR INTEGRATION 


1. Swmmation is the process of finding the sum of any number 
of terms of a given series. This can be accomplished either if the 
law of the series is known or if a sufficient number of terms is given 
to enable the law to be ascertained. As will be shown in Chapter XX, 
if no mathematical law is apparent, methods can be applied by 
which the approximate sum of a series can be obtained. 


2. Consider a function F’(#) whose Ist difference is f(x). Then 


we have 
FQ) -F(0) =f() 
F(2) -FQ) =f() 


FQ@) —F@-1)=f(a-1) 
Pa+))-F(@) =f) 


EF(n-1)-F(n—2)=f(n—2) 
Fin) —F(n—1)=f(n—1). 
Summing both sides, we obtain 
F(in)-FO)=f(0)+fQ) +...+f(n—1) 
or F(a)-F(O)=f/(0)+fQ)  +...+f(a—1) 
or Fn) —- F(a) =f(a)+fat+1)+...+f(m—1) ...(). 
It is clear, therefore, that the sum of any number of terms of a 
series of values of f() can be represented by the difference between 
two values of another function F(a) whose Ist difference is f(a). 
By analogy with the system of notation already adopted for ex- 
pressing orders of differences, the process of finding the function 
whose Ist difference is f(x) may be denoted by A f(a). It is 


customary to express f(0)+/(1) +... +f(n—1) as = J (@), the 


terms at which the summation is commenced and terminated 
(designated respectively the inferior and superior limits of sum- 
mation) being indicated in the manner shown. 
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3. The process of finding the value of #’(«) is known as Finite 
Integration and F(x) is called the Finite Integral of f(x). Where 
the limits of summation are known we obtain by summation of 
Ff («) the Definite Integral of f(x); 1f the limits of summation are 
not expressed we obtain merely the Indefinite Integral of f(z). 


4. As stated above, in obtaining the indefinite integral of f(«) 
no point is specified at which the summation is to commence and, 
since an unknown number of terms of the series is included, it is 
necessary to include in the value of #’(«) a constant term which is 
of unknown value. 

This constant vanishes in the case of definite integrals, since if 


Sf (0) = F(a) +e, 
shee = f(a) =[F(n) +0] -[F (0) +e] 
= F(n) — F' (0). 


5. It is obviously always possible to find the first difference of 
any function, but it does not follow that every function can be 
integrated. The functions which can be integrated are limited in 
number and the process of integration rests largely on the ingenuity 
of the solver aided by such analogous forms as may be obtained by 
the formulas of finite differences. 


Thus we have Aa* =(a— 1) a’, 


whence it is easily seen that 


Aa* 
aes 
oo alg 
Ea 

and therefore, since the result of differencing _ i is to give a%, 
we have 

a” 

Yat = +, 
a—1 


where c is the constant introduced by integration. 
The sum of the series a’ + a’? +... + a’+”— is at once obtained 


r-+n-1 
by finding the value of the definite integral 2% a*, which by §2 


rT+n 


, , 
is equal to <— -—, which agrees with the familiar result for 


the sum of a geometrical progression. 
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6. Similarly since Aa™ = mar), 


(m-+1) 
by analogy Lal) = = ro ad 


7. The following table gives the values of some of the simpler 
integrals. They should be verified as an exercise by the student. 


Function Indefinite Integral 


4 ae ee 
a* 
ae oa 
a (m+) ; 
— m+1 
al-™) alee, 
= (1) 
— a an +b) (ax—1+b)...(av—m+b 
(az +b) (ax —1+ 6)...(az—m+1+6) See Nn tte 
a eta as (ax—1+46)(ar+b)...(ax+m—1+b 
(ax +b)(ax+1+45)...(ax+m—1+45b) a a 
1 1 
(ax + b)(az+-1+b)...(az+m—1+5) —a(m—1)(ax-+b)(ac+1+b)...(az¢m—2+b) °° 


i I 
(aa +b)(az—1+b)...(az—m+1+b)) -—a(m—1)(az—1+4+5)... ai 


8. If the form of the function is unknown, a general formula for 
the sum of a series of values may be obtained as follows, since 


f(a) = (0) + eA f(0) +22 wf) 
eae 

Integrating both sides, we have 

Sf (e)=0-+0f(0) +2 2—V a0) + 22-DE— ayroy+..., 

or, integrating between limits, we have 


Sy@-= nf (0) pS DA yo) s Se) aero) 4 a 
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or, more generally, 
“ES f@)=f()+ f+) +... +fl@en-2) 


= nf (a) + @ 2D pp(ay PO“ DO Darra) 4... 


It is instructive to obtain this result by the method of separation 


of symbols. For 
F(a) =f (a), 


Fai) =f ORO) 
flat 2)=(1+ APS (a) 
Kien net ayy 

Adding 


S f(a) ={1+(14 A) + (14 AP 4... (4A) f(a) 


_(i+A)" —1 
~ (1+A)-— @¢ayci’@ ») 
={14nag "are... 1} 


- f@) 


nn a OR are LF @ 


=nf (a+ Vas (SOS ee 


9. An example will illustrate the use of the formula. 


Example 1. Find the sum of the first n terms of the series whose 
initial terms are 1, 8, 27, 64, 125. 
We have 
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Whence 
S f@=nx142G— x a el 12 
oe ee 


_n(n+2n?+n)_ ni(n+1) 

. 4 ~ ao” 
which agrees with the formula for the sum of the cubes of the 
natural numbers. 


10. Where it is desired to integrate a function which is the 
product of two factors, the following device may often be utilised 
with advantage. 


Let the function be Y = UzgVz. 
Then Ag Vz = Ug Vag. — Uxe 
= Ups (reg — Uz) + Ug (Uz+1 — Ux) 
| = Uggs, AVz + Vz Atz. 
Integrating both sides of the equation, we obtain 
Uy Ve = LUzy, AVy + Dz Auy 
or Bivig Atte Cet Bite - loess ects... (8). 


Thus, if the original function can be put in the form v,Auz, its 
integral can be made to depend upon that of uz,,Av,, and, if the 
latter is in a form which can be readily integrated, the value of an 
apparently intractable integral may often be obtained in this way. 


Example 2. To find the value of 2a’. 


‘ a* 
Since 2a” = ———,, we may write 


a—]’ 
Hy) 
Lca® = La ae , 
a-l 
Using the above formula (3) we get 
Aa® aa* qt 
20 ani gai 
e+. 
= an 3 = , since Av = 1, 
a-—l a-—l 
xa” an 
“s-1 @oep * 


H. 1. B. 1. 4 
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11. Sometimes it may be necessary to apply the formula more 
than once in order to reduce the integral by stages to a standard 
form. The process is illustrated in the following example. 


Example 3. To find the value of 22*2%. 
Remembering that A2* = 2*, we may write 
22% a8 = Las AQ? 
= 92% — 22711 Ag? 
= 92% — S241 (B22 + Ba +1). 
It will be observed that in applying formula (3), the degree of 
x, in the terms within the integral, has been reduced by unity. 
Proceeding as before we obtain 
22708 = a8 2% — 3 (3a? + 8a 4+ 1) A2*+) 
= 982% — (27+) (82? + 8a +1) — 227°A (32° + 824 1)] 
= 2" (a — 6a? — 6a — 2) + 22°** (6x + 6) 
= 2% (a — 62° — 62 — 2) + > (62 + 6) A2*+# 
= 2° (x — Ga? — 6a — 2) + [27+* (6x + 6) — 227A (6x + 6)] 
= 2° (a* — 62° + 18x + 22) — 227+? x 6 
= 2° (a? — 6224+ 184 + 22)-—6 x 27 +6 
= 27 (a? — 6a? + 182 — 26) +c. 
The above process is analogous to that of “Integration by Parts,” 
which is dealt with in the Integral Calculus, Chapter XVIII, § 6. 


CHAPTER VIII 
FINITE DIFFERENCES. DIVIDED DIFFERENCES 


1. A simple method of interpolation is available, where the in- 
tervals between the given terms are unequal, by the method of 
Dinded Differences. 


2. The application of the method rests upon the assumption, 
which, as has been shown, is the basis of all theorems for interpola- 
tion by means of Finite Differences, that f(x) 1s a rational integral 
function of x of the nth degree. 

On this assumption, it can be shown that f(x) can be expressed 
in the form 

A,+ A,(@ — a,) + A,(a—a,)(@ — a4) +... 

+ A, (a — a,)(@— ag)... (L&— An), 
where A,, A,,... An, G, Ge, ... Gy are constants. 


3. In order to apply this formula in practice, it is convenient to 
introduce a scheme of notation on the following lines, where the 
symbol of operation is denoted by A’ in order to distinguish it from 
the ordinary A. 


F()-F(0) | A’ f(ar)—A'F(0)| 42 F (ay) - 4/2 f(0) 
ay 


f(a2)-F (ca) | A’ f (a2)— Af (a)| Af (a) - Af (mn) 


aaa 


ag— a, * a3—a, 


f(as)—F (42) | AF (as) — A’ f (a2) 


ag— a2 a4— Ag 


Las) — F (4s) 


@—A3 


Generally, we have 


™ A”™ f(a, —A’-1 ™ 
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Hence (a) =f(0)+ a4 f(0), 
and SF (a2) =f (th) + (2 — 1h) A’ f(a) 
=f (0) +a, A F(0) + (a2 — a1) {A’F(0) + a, A? f (0)} 
=f (0)+a,A’f(0) + a, (a,— a,) A” f(0). 
By proceeding similarly for further terms, we find that we can 
write generally 
f (0) =f (0) + 2d’ f (0)+2(@—a) A? f(0) 
+a (e—a)(@— a) A*F(0) +... (1). 
This general form can be readily established by the method of 
induction. 
By giving appropriate values to a,, d:,... the ordinary formulas 
applicable to equal intervals can be at once deduced. 


4. The genera] method of working will be shown more simply 
by an example. 

Example. Find the value of log 4°0180, having given the fol- 
lowing data: 


Number Logarithm 


*6020600 
6034367 
°6045824 
°6052835 
“6060480 


Transposing the origin, we have the following scheme: 


x f (2) A’ A’? 


0013767 — 000317 _ 


00000 | 6020600 Dig = 108402 | Sas = - 01186 
00127 | -6034367 ett 108085 , =~ 0130 

| 0-0233 | -6045824 rae = 107862 ae 70187 
0:0298 | -6052835 aan ‘107676 


0°0369 “6060480 
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We have to find f(‘0180), which is, by the above formula, 
F (0) +0180’ f(0) + 0180 x °0053A7 (0) 
[the further terms will not affect the seventh place of decimals], 
where /(0) =°6020600, A’ /(0) ="108402, A” /(0) = — ‘0136. 
Thus log 4°0180 = -6020600 + 00195124 — ‘00000130 = 6040099 
to seven decimal places, which agrees exactly with the true result. 


CHAPTER IX 


FINITE DIFFERENCES. FUNCTIONS OF TWO VARIABLES 


1. Questions involving functions of two variables arise frequently 
in actuarial practice. Thus the tabulated values of functions (e.g. 
annuities) dependent upon two lives may be given only for com- 
binations of quinquennial ages in order to economise space. If the 
value corresponding to any other combination of ages is required, 
resort must be had to methods of interpolation. 


2. In considering the problem of the changes induced in the 
value of f(z, y) by finite changes in the values of # and y we must 
consider z and y as being independent of each other. Clearly, if 
y were a function of x the expression f(a, y) could be made to 
assume the form of a function of « alone and it could be dealt 
with by the methods already developed in previous chapters. 

Thus x may vary while y remains constant, so that, if « changes 
to +h, the value of the function becomes f(a +h, y); or y can 
vary while # remains constant, giving a value f(x, y +); or both 
a and y can vary independently, giving a value for the function of 
S(ath, y+). 


3. We shall proceed first to discuss the problem where the 
values of the function are given for combinations of successive 
equidistant values of # and y. 

Thus we may have 


F@Y) — fethy) flar%hy) — ...f(e@+mh,y) 
f(aytkh) f(ethyth) f(c+2h,yth). ...f(a+mh,y +h) 


F(a,y+nk) f(ath, y+nk) fa haat) we f (a+ mh, y+ nk) 


As has already been seen in the case of functions of one variable 
(Chapter ITT, § 6), this scheme can be simplified, for the origin can 
be placed at the point (a, y), and the unit of measurement can be 
taken as h in the case of the variable « and & in the case of the 
variable y. 
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The scheme then becomes 


F(0, 0) fA,.D PEO AM lees sc tas F(m, 0) 
F(O, 1) FQ, 1) C0 ee eee F(m, 1) 
F(0,n) fl, n) if Came) Bites 2, F(m, n) 


4. Since « and y may vary independently, a fresh scheme of 
notation must be introduced to express the variations which may 
arise. Thus A, will be used to denote the operation of differencing 
with respect to #, y remaining constant, a corresponding significance - 
attaching to A,, so that 


A, f(0, 0) =f (1, 0) = A0; 0), 
A, f (0, 0) =f (0, fy —F (0, 0), 


or, using the method of separation of symbols, 
F(1, 0) = (+ Az) f(0, 0). 
Accordingly, we have | 
Sf (m, n) = (1 + Az)™ (1 + Ay) f(0, 0) 
=(L+mAz+($) A+...) (1+ndy + (5) d%+...) £©, 0) 
= f (0,0) +mAzf(0,0)+ (3) A% £0, 0) + (3) A%.£(0, 0) +... 
+nAyf (0,0) +mn AeA, f(0,0)+ (5) nA%Ay f(0,0)+... 
+(3)4% £(0, 0) +m (3) ArA%,f(0,0) +... 
? + (5) SFO, 0) +... 
Here A?,, A3,, ... - be written down by differencing the rows of 


the table of the function; similarly A?,, A’,,... are the differences 
of the columns of the table. 
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To find A,A,, A?,A,, etc. we have 


4, A,f(0, 0)= Az [f(, 1) —f(, 0)] 
=f(1,1)-f (0, 1) -f(1, 0) + (0, 9), 
A?,A,f (0, 0) = A*%[f(0, 1) — f (0, 9)] 
=f (2, 1) — 271, I). +f (0, 1) —f(2, 0) 
+ 2f(1,0) —f (0, 0), 


and so on. 
Example 1. Table XVI of the “Short Collection of Actuarial 
Tables.” To find A,},,,, having given 
eh = L1LG6D, Ad on = 13190, A gh igo = °15494, 
Ag .x, = 09809, Ag) eg = °11089, 
Ag my = 07812. 
Here m=}, n= 3, and 
By, » 1521, At,= ‘00783, A, A, = — ‘00291, 
A, = — ‘01860, A?, = — ‘00137. 
Whence A351, = ‘10776, the correct value being ‘10773. 


5. An obvious method of procedure involving only first differences 
is as follows. Obtain the value of f(0,n) by interpolation between 
the values of f(0,0) and f(0, 1). Similarly, obtain the value of 
J (1, n) from the values of f(1,0) and f(1,1). Finally find f(m, n) 
by interpolation between f(0, n) and f(1,7n). Thus 

FO, n)=1—nf(0, 0)+ nf 0, 1), 

FQ, n)=1—nf (i, 0) +f, 1), 

Ft (m, n)=1—mf (0, n)+ mf(1, n) 

= f (0, 0)+ mA, f (0, 0) +nAy f (0, 0) + mnd,A,f(0, 0) 


=I—-mi_-nf(0,0)+n.1—mf (0, 1) 
+m.1—nf (1, 0)+mnf (1, 1) ...(8). 
Employing this formula in the example given above we find 
stm = 10822. 


The method is suitable if only a rough approximation is required, 
but cannot be depended upon to give an accurate value. 
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6. Obviously the method can be extended by taking higher 
orders of differences. The disadvantage of this procedure is that 
it involves the calculation of further values of the function corre- 
sponding to a given value of x as a preliminary to applying the 
interpolation formula to find the value of f(z, y). The process 
thus becomes laborious and moreover we do not necessarily obtain 
identical values for f(a, y) if we interpolate first with regard to « 
and then with regard to y, or vice versa. 


7. As in the case of functions of one variable, we shall expect to 
obtain the best results when the principles of central differences 
are applied, 1.e. when the required term occupies as ‘nearly as pos- 
sible a central position among the terms employed in the formula. 
The difficulty is that, in dealing with functions of two variables, 
we cannot adapt our formulas to any system of values which may 
be given. Thus an inspection of the advancing difference formula 
(1) shows that it involves points whose coordinates form a triangu- 
lar plan which may be illustrated thus: 


0, 2) 
(0,1) G1) 


(0,0) (1, 0) (2,0) 


This illustrates the formula where two orders of differences are 
taken into account, the black dot representing the interpolated 
term. It will be seen that the scheme is hardly satisfactory from 
the point of view of central differences. For most practical pur- 
poses, however, where ordinary actuarial functions are involved, 
formula (1) will give satisfactory results. 


8. Formulas embodying the principles of central differences can 
conveniently be obtained by an adaptation of Lagrange’s for- 
mula. This formula applied to functions of two variables has not 
the same wide application as the ordinary formula of Lagrange 
previously given in Chapter IV, but, as will be seen below, it gives 
expressions for f(x, y) in terms of the neighbouring values. 
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9. General formula for 4 points. 
Taking all combinations of two terms except those which give 


rise to 2 and y?, let 
I (@, y) =A (2@—B)(y—b) + B(x — B) (y — a) + C (x@— a) (y—a) 
+ D(a— a)(y—6), 
then I (a a) =A (a— B)(a—5), 
Sf (a, 6) = B(a— 8) (6 — a), 
(8, a) =D (8 —a) (a—6), 
f (8,b) = C(B—a)(b —a), 


whence, substituting for A, B, C and D in the original formula, 


2 — =) 
SF (@, y) =f (a 0) + FO, ee ae 


(x —a)(y—a) a) (y — 6) 
+f (8, ) Baayboa) 17 A 0) eae 


10. General formula for 6 points. 
Taking all combinations of two terms, let 
f (@, y) =A (@— B)(y — 6) + Bie B)(y— 4) + C(w—4@) (y—a) 
+D(w—a)(y—b)+ E(w —a) (@— 8) + F(y—a)(y—3). 
Taking the points a:a,a:b,a:c, 8: a, B:b and y:a, and pro- 
ceeding as before, we arrive at the result 
—B)(y—b) | (@-a)(@—B) , Y¥-a(y-4) 
He ~Fe0) (gay baa)* @=w= a) * OaaNera) 
z—B)(y—a) , (y — a) (y ~ 6) 
- flab) {e-EN —a)(b— Da) $e 
+ f (4,0) oe =) a — 5 
J a” a)(y—b) | (w—a)(a—8) 
AB0)\G—ayo=a) * B=HG=AN 
+fBD gg 9 
on — 4) (#— B) 
+249 ay y=B) 
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11. General formula for 9 points. 


Taking all combinations of two terms, each involving 2, with 
two terms each involving y, let 


J (@,y) =A (@— B)(z—)(y—b) (y—c) 
+ B(a— B)(@—-)(y— 4a) (y—¢) + C(@— B) (@—) (y — a) (y—b) 
+ D(a — a) (a —)(y— 6) (y —¢)+ BE (@—a)(@—-y)(y—a)(y—c) 
+ F(a — a) (@— y)(y— 4a) (y — 6) +G (@—a) (2 — B)(y —b) (y —0) 
+ H (a —a)(@—8)(y—a)(y—c) +I (@—a)(«— B)(y—a)(y—}). 


Whence, proceeding as before, we have 


- (x — B)(@—) (y—6) (y—c) 
POD=LOD @=Bya=y)(a=B)(a—e) 


(x — B)(a—-)(y —a)(y —c) 
THO C= ay@— yO Salo 


(x — 8) (w—y)(y —a)(y ~8) 
+ H(t °) = By(a—y) (0a) (08) 


(@—0)(@—y)(y—B\y—o) 
tO Ds (C= aie OOO) 


(7@—a)(e—Y)(y— a) (y—C) 
Pa aaa =ag( aa oeo 


ox Demag akg 8) 
tI Gl am Gae 9) =a) 


(a —a)(x— 8) (y—b) (y—c) 
MeCN = (a (CET 


(x —a)(a— 8)(y—a)(y—c) 
tT On eal aa (0) 


(w—a)(x— B)(y—a)(y—5) 
+ f(y, Cio Ga a eieidais'e ca cian (6). 


12. Formula (4) is the general formula corresponding to the 
method of § 5; by altering the notation the identity of the two 
formulas (8) and (4) is apparent. 

Formula (5) includes six values of the function, the co-ordinates 
being related in the manner shown. It will be seen that the formula 
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can be applied to any of the following groups of values, the black 
dot representing the interpolated value: 


(0, 2) (0,1) (1,1) (0, 1) 
fo) (o] [e) (o) 
(1) (4,1) (-1,0) (0,0) (1,0) (~1,0) (0,0) (1,0) 
(o) (o) (eo) fe) (e) fe) (eo) (o] 
(0,0) (1,0) (2,0) (0, 1) (<1; = G,-1) 
(o) (eo) (o) (eo) (e) [o) 
(1) (11) (aii) 
(1) (4,1) (1, 1) 
(-1,0) (0,0) (1,0) (-1,0) (0,0) (1,0) 
fe] (eo) (eo) fe] (e] (°) 
ae) -y-> ep 
(iv) (v) 


Other systems could be written down, but the above are suf- 
ficient for purposes of illustration. 

System (i) is obviously the same as the advancing difference 
formula. It is obtained by writing a=0, 8=1, y=2, a=0, b=1, 
e=2 in formula (5). 


Example 2. Find A,}.., from the data of Example 1, using 
the Lagrange formula applicable to system (i). It will be found 
that the same result is obtained. 

System (i1) should be expected to give a formula which will be 
more accurate than the advancing difference formula, since the 
interpolated term will occupy a more central position in relation 
to the terms employed. The formula is obtained by putting a= 0, 
B=1, y=-1, a=0, b=1, c=—1 in formula (5). Working on 
the same example as before and taking the origin at the point 
(30, 60), the values of the function entering into the formula are 


Ag}. =10080, Act.gs='09809, Agi .go = ‘11669, 
30:65 = "14006, <Asi.65='11039, A3}. 69 = 18190. 


We obtain as a result A,},,..=°10770. The degree of approxi- 
mation, though close, has not, in this instance, been improved. 
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System (111) is an inversion of system (ii) and should be useful 
for interpolation where x and y have negative values. 

The lack of symmetry of systems (iv) and (v) suggests that they 
are not likely to yield good results in practice. 


13. When nine points are used, as in formula (6), the system is 
represented by the following diagram : 


(=1, 1) (0, 1) (1, 1) 
(—1, 0) (0, O)e (1, 0) 
ce) ° ° 
(Ah |) (0, =) (1, =) 
° fe) (°) 


It will be seen that this scheme embodies all the principles of 
central differences and should therefore give good results. 

Taking the previous example with the origin at the point (30, 60) 
the six values entering into the formula for system (ii) are used 
together with the following additional values : 

Adirgs= 08435,  Aghgg = "12132, Ab gg = “15972. 

Making use of formula (6) the interpolated value is found to be 
‘10771, a slightly better approximation to the true value than 
those obtained previously. 

On general reasoning we should expect a somewhat better result 
by taking the origin at the point (30, 65) so that the interpolated 
value would occupy a more central position. The values A,j.5, 
Asi :55 Agk.g5 entering into the immediately preceding calculation 
are excluded, and the following values introduced : 


As 279 = 06642, si: = 07812, 3h 79 = (08756. 
On working out the result, however, we arrive at the value 10848, 


which is a worse approximation than the value obtained by the 
rough method of § 5. 


14. This apparent inconsistency illustrates one of the chief 
difficulties of interpolating between functions of two variables, 
namely, that one does not necessarily obtain a better degree of 
approximation by proceeding to a higher order of differences or by 
employing more terms in a formula. Changes in the value of f(a, y) 
occasioned by alterations in the values of # and y may be so con- 
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siderable that distant terms may have such a disturbing effect 
upon the formula used as to upset the agreement between the 
approximate interpolation surface and the true surface which 
represents f (2, ). 

It is thus difficult to say what will be the degree of approxima- 
tion of a given formula, but an inspection of the course of the 
differences will be some guide as to the advisability of introducing 
further terms into the calculation. 


15. Other devices may sometimes be adopted which enable the 
interpolation to be reduced to the work of a single variable inter- 
polation. 

Thus, if the sum of a and y is a multiple of 5, by suitably 
selecting the origin we may write 


Sf (a, — @) =f (0, 0)+a[ fd, ail 1) —f (0, 0)] 
+ (3) (F@ -2)- 20, -)4F@, 0] 


at — eg —% 


3 0)-—a#.2—2f(1, —1) 
+25" F (2, =2) ~ (7). 


By referring to the point diagrams on previous pages it will be 
seen that the process is equivalent to interpolating along a diagonal 
line running through the various points. The formula is of the 
advancing difference type; the corresponding central difference 
formula would preferably be employed in practice. 

J. Spencer has given (J.J.A. Vol. 40, pp. 296-301) examples of 


the use of several ingenious methods of this character. 


CHAPTER X 


DIFFERENTIAL CALCULUS. ELEMENTARY CON- 
CEPTIONS AND DEFINITIONS 


1. In the subject of Finite Differences we were concerned with 
the changes in the value ofa function consequent upon finite changes 
in the value of the independent variable. In the Differential Calculus 
we consider the relation of Ay to Az when the value of Ax is made 
indefinitely small. 

The application of the Differential Calculus is largely limited to 
such values of a function as are finite and continuous, and, unless 
otherwise stated, this limitation is to be implied in the following 
demonstrations. In practice these conditions are almost universally 
fulfilled by functions entering into actuarial calculations. 


2. Let y=/f(#) and let « receive an increment h. Then the 
change in the value of y is measured by f(x+h)-— f(x) and the 


rate of change of y es . The limit of this expression 
when h->-0 is called the Differential Coefficient or First Derived 
Function of f (a) with respect to a. 

The operation of obtaining this limit is called differentiating f(a). 


Using the notation of Finite Differences the differential co- 
efficient becomes 


= 
ser oAe 
and is variously denoted by “u f’ (a), f (2), — Die). 
The symbol = or its equivalent represents an operation of the 
character described; the elements dy 
and dx must not be regarded as separate y Q 
small quantities. P 


R 


3. The geometrical representation of 

the differential coefficient is illustrated 

in the accompanying diagram. © fF MN X 
Let the curve shown represent the 
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equation y= f(x). Let OM=x2 and ON=a+h, and let PM and 
QN be the corresponding ordinates. Let PR be the perpendicular 
from P on QN and let QP be produced to cut OX at T. 


Then Oe ae eee 


When the point Q moves up to, and ultimately coincides with the 
point P, the line QPT becomes the tangent to the curve at the 
point P. The limiting value of a! is therefore the 


tangent of. the angle which the tangent to the curve at the point 
(z, y) makes with the axis of a. 


CHAPTER XI 


DIFFERENTIAL CALCULUS. STANDARD FORMS. 
PARTIAL DIFFERENTIATION 


1. The differential coefficient of any particular function can, of 
course, be obtained by direct calculation, but the process can 
usually be simplified by the application of the following general 
rules. The general similarity to the propositions already demon- 
strated for Finite Differences will be apparent. 


I. The differential coéfficient of any constant term is zero. 
This is evident since a constant is a quantity which does not 
change in value in any mathematical operation. 


II. The differential coefficient of the.product of a constant and 
a function of « 1s equal to the product of the constant and of the 
differential coefficient of the function. 


Thus =F fe.f@= Lr PEF“ F 
ip fet®)-f@ 
h->0 h 
df («) 
we ee (1) 


III. The differential coefficient of the algebraic sum of a number 
of functions of x is the sum of the differential coefficients of the 
several functions. 


Let y=u+v+w-+..., where u, v, w, ... are functions of x, 


then Ay = Au+ Av+ Aw+ ... 
Ay Au, Av , Aw 
and ; in fo ae 


which, by proceeding to the limit, becomes 
dy du dv dw 


ae as “ae ae Coes rer essere nsees (2). 


IV. The differential coefficient of the product of two functions ts 
the sum of the products of each function and the differential coefficient 
of the other. 


H. T.B.I. 4) 
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Let y = ut, 
where wu and v are both functions of x. 
Then Ay =(u + Au) (v + Av) — uv 
= udAv+vAu+ Audyv 
= uAv + (v + Av) Au 
Ay 
and AY a SY 4 (v + Av) SY 
whence, taking the limit, when v + Av, 
dy__ dv, du 
ae =U ae +v Zz RC OC Cee (3), 


which may be written 
ldy idu_ ldv 
ie SS ee (4). 
This result may be extended to include the product of any number 
of functions. 
For if y= www; let vw = z, then y = uz. 
ldy_ldu ldz 


Whence y + = w ae ae 3 ax . 


But ldz_ 1ldy 1 dw 
zdz vax wdx 
ldy lduldv 1 dw 


Therefore pao BE +o atwds cc (5). 


Multiplying by wow, we obtain 
OY a ww Se + wn 2 + ww stoi WS vowel (6). 


Similarly for the product of any number of functions. 


V. The differential coefficient of the quotient of two functions 1s 
(Diff. Coeff. of Numer.) (Denr.) — (Diff. Coeff. of Denr.) (Numr.) 


’ Square of Denominator 


Let =, 
: v 
" u+Au ui vdu-—udAv 
i a eae v(v + Av) 
Au . Av 
v— —u— 


Ac v@+ar) ’ 
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whence, taking the limit, 


which may be written 


yda uda vda 


VI. The differential coefficient of y with respect to x, where y is 
a function of u and u is a function of x, 1s the product of the dif- 
ferential coefficients of y with respect to u and wu with respect to «. 


Ay Ay Au 
For As = Ad . a ’ 
whence, taking the limit, 
dy _ dy du 
da du'dg Us (9). 
oe dy dy du dv - 
Similarly dc du'do'de cna (10), 


and so for any number of functions. 


2. Various standard forms can now be developed, mainly from 
first principles. It is instructive to note the points of analogy with 
the corresponding forms for Finite Differences. 


(3) = a”, 
dy _ = (a+ hyn — a” 
da h—>0 h » 
lett 
= Lt : 
h->0 y h 


Expanding by the Binomial we have 


eet at Get ~| 


= Ltn E ae hy 
h->0 2 £ 


3 


=n), | 
5—2 
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2. 
=a Lt i [2 + hlog, a + (log, ay+... -1| 
h 


=a’ Lt log, a+ : (log, a)? +. ‘i 
h->0 


= a” log. a. 
aig GS a 
If ye, ld. lead 
(111) y = logaa. 
Then ay = x. 
d(a¥)_ d(a¥) dy 
But ap = dy : ide 6 


Hence, using the result established in (ii) and remembering that, 
aY =a, we have 


ax dy 
dag OY Bea a 
or it = wlog.a. SE. 
Whence ees. 
dx sxlog.a 
oe oe 1 a 
Tr y = log. a, da s«logee x’ 


(iv) y=([f(a); y=e™; y=log. f(z). 
If y¥=([f(@)), 
dy_a[ f(a) df(a) 
dz  df(a) ° dx 
=n[ f(a) .f' (2). 


Similarly if y=e%, Wes. ¢(@), 
; d 4 
ove y=log f(a), = on 
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(v) y= sin 2, 
dy _ sin(a+h)—sin x 
as 
2 sin 5 cos (1 + *) 
= Lt oe ae 
h-0 
sin - 
2 h 
+e COS (= + 5) 
p. 
= COS &. 
(vi) y =sin— 2. 
Then sin y = &, 
whence Fyn cs y= V1 = sinty = V1 a 
ay 1 1 
erefore es aaa 
dy 


3. The values of differential coefficients for the other trigono- 
metrical functions can be found by methods similar to those employed 
im (v) and (vi). The results are given in the table below and should 
be verified as an exercise by the student. 


. Differential 4 Differential 
Function Coefficient Function Coefficient 


a nmagn—t 
a a* log, a 
e* e* 
aod 
a log, a 


sin~!4 


cos~!z 


log, £ 
tan—! 4 
log, # 

; cot~lz 
sin x COS £ 


COS & — si 
Be sec! a 


tan £ sec? x 

= 2 
cot z cosec* £ cosec™! ¢ 
sec £ secxv.tanz 


cosec £ —cosec #.cotz 


70 DIFFERENTIAL CALCULUS 


4. Logarithmic Differentiation. This method is of special value 
in two cases. Thus if y= uvw..., where u, v, w,... are functions 
of z, then ; 

log y = log u + logy + logw+..., 


ldy ldu ldv, 1 dw 


— y deh de. ve a det 
dy _ - ldu ldv 1 dw 
i da luda vada’ wdz TE: 


a result which agrees with that already obtained in § 1. 
Secondly if y= u*, wu and v both being functions of «, 


log y = v log u, 
ldy_vdu dy 
and in eae * un 


or 2 
1 


5. We will now give some saisooMatrents examples of differen- 
tiation. 


: a+ x 
ie acl e 
By the ordinary rule for a quotient 
\ ») 2 (a +2") 2 2(6 + 2) 
+) 01 


dy _ : 
le (64a) 
Sa) (a+ a*) a+ 2ba-a 
| (b+ay Cray 
a _* a Meal .| | 
y ye 1 — 5a + 6a*” 


This can best be treated by resolving the expression into partial 
fractions. Then 
« 4 8 : © 
a ae ey | 
dy _44(1 — 3x) a (I — 32) _ dL = 20) ~ 2a)" d(1 — 2a) 
dz ~d(l—32)"° da d(1—2x)° da 
12 10 
~ (1-32) (i - (1 — 2a)?" 


: 
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(111) y= cee 


(iv) 


: a. 1 af 
Ae 4)  —_ 

(tan Ve Ve—1 d(s*—1) 
‘d(@—1) dx 


= ——.- 4-17 ba 


“ ava? —1- 


(vi) y= (= 2 


1, —~= 1 
log y = = log a — 2 —— log a — 8, 
1 


1dy _1d (log #— 2) 2) d(a—2) ae 
y da — a? a areerss es 
1 


— 


log #—-3.-=, 


ll 1 a ge 2 
Be — 2): ae 2) ia ba 8. | 
: i 2th 1, x2—2 , 

| 


dy” {[“2—2 
Whence dy __ (ex 3; Ee = 8) aD a! SF 
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(vii) y= ar. 
log y = x log a + c* log b, 


ldy _ 
ge Oe ey, 


oy a* .b® (log a + c* log c log b). 
1 
(viii) y =a + oF, 


In this case logarithmic differentiation must be used, but for this 
purpose the two terms must be taken separately. 


1 
Let at =u and a=», 
Then dy du. dv 


da da. dx’ 


log u= « log a, 


du 


dg & (1 + log 2). 


Also log v= : log a, 


dy er) 
Therefore ine (1+ log 2) +a* “(1 — log 2). 


(ix) Differentiate log, x with regard to a”. 


Let y=log,« and z= a’. 
dy dy dx dy 1 
Then dz da'dz~ da’ da’ 
dx 
Therefore dy Lit 
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6. In dealing with cases where a function of two variables is 
involved it is convenient to adopt methods similar to those used 
in Finite Differences (see Chapter IX, § 4). Thus we define Partial 
Differentiation as the process of differentiating a function of several 
vaniables with reference to any one of them, treating the other 
variables as constants. 


This process is denoted by the symbols ke =, etc. We will also 


ox’ 
use the symbol 6x to denote a small change in the value of «. 
Let u = f(a, y). 
Then ut+du=f(«e+h, yt+k) 
and 


bu=f(e+h, y+h)— f(a, y) 
_f@thy+h-f@yth ,  f@yth-f@y , 
; | : LB. 


Proceeding to the limit when h and & successively — 0, 
| th,y+h)-f(eyth_ 


=< a (x, y), since k > 0. 


Ou ( 
Therefore ou = aa bx + a by, 
du du du by 
and =. = Ae + ay a “ 
du du , ou dy 
or, when 6x0, Pte hl oe (11). 
tT. If f(e, y) =9, 7 = 0 and 
au , du dy 
On Oy da 
Ou 
dy Ox 
eee 12). 
Whence ae —o (12) 
oy 
Example. If a+ cy+y' = 0, 
ou ou 
then a 20 + Y, een 
dy +, 2e+y 


and 


dx «2+2y 


CHAPTER XII 


DIFFERENTIAL CALCULUS. SUCCESSIVE 
DIFFERENTIATION 


1. When = is differentiated with respect to 2, we obtain a 


further function of « which is known as the second differential 
coefficient or second derived function of y. This operation represented 

d /dy\ . d?y ” - gg 
by aa (G4) is generally written Tat OF f ip). If the ae 1S 
differentiated n times with respect to 2, the result is called the nth 


e e e = e Ld : = d” 
differential coefficient or nth derived function and is written — 


Other symbols for the nth differential coefficient are FO (&), DY; Un. 


2. In many cases the value of the nth differential coefficient can 
be found readily by inductive reasoning. 


Eaample 1. 
y = log (a+). 
| ag eee 2 J 
oS py Oe pasa? TS (a +a) ; 
and by analogy . 
Ge 1)! 
n= — | ye! , 
Yn =(—1) @ray 
Example 2. ,, 7% 
y=atbat+ ca®+ datt+...+ka, 
W= b4+2cat+ 3da*+... +n. kar, 
Ye = 4 2c+2.3dax +...¢n(n—1) ka", 
Yn = y nik. 
Example 3. | 
Taz—1 


IT obepeey eM 
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Making use of the method of partial fractions, 
4 5 
= 73a ~ Toe 
4.3 a. 2 
A> ion TB) 


5 4.3%.n! 5. 2° mt 
eBay (12a 


3. Leibnitz’s Theorem. 
This theorem gives an expression for = : where y is the product 
of two functions of z, say wu and v. 

‘It should be noted first that, as shown in Chapter XI, § 1, the 
operator - or D obeys the distributive. and index laws and is 
commutative in regard to constants. In these respects it is similar 
to A (see Chapter II, ya 


Now , £ wav tu, 
and we may therefore write — 
D (uv) = D, (wr) + D, (ur) 
= (D, + D,) uv, 
where D, operates only on wu and differential coefficients of «w and 
D, operates only on v and differential coefficients of v. 
Therefore 


D* (uv) =(D, +D," uv 
=| Dr +nD 2. LD a(t ) Dr Det .. + Dy] we. 


Now 
tay 
Le ate "a 
a 1 (t?\ dv Os 
DoDy (ue) = Di (uo) = Faia 
ete. 
Hence 
ad™(uv)_  d®u dv d™ tu | (n\ du d™*u d"y 
a oontn da’ da (3) dat’ dara to + Gan 
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This formula can also be established by induction, but the proof 
by this method is left as an exercise to the student. 


4, Examples of the application of the formula are given below. 


Example 1. y = xe, 


i Oe 
Since ies 
take 2? as the factor v, and e® as the factor wu. 


Then gh = 32°, etc. 


da 
Therefore 
Yn = Care +n. 320%. a" em + (5) 6a. a®—2 07 4. (") 6. a™%e%, 


can be written at once as a”e™, it 1s convenient to 


Example 2. If y = x" (log x)’, 


prove that 
LY nse + LY n+1 = 2 x n!. 
We have Y= an”, aone + (log x) n. gr 
or xy, = 2x2" log x + nx" (log x) 
aa a ee ee ee (1). 


Differentiating both sides of the equation 
Dyn 
LYo + Y; = — +2 log a.nx"™ + ny, 


or ay, + xy, = 2x" + 2nz" log x + nay. 
Substituting for x" log # from equation (i) we obtain 
Ly, + LY, = 22" + nxzy, — ny + nxy, 
or Buy — (Qn — 1) ay, + Py =QA™ ......ecreeeee (11). 
By differentiating each term of this equation » times, making 
use of Leibnitz’s Theorem, we arrive at the result 
Unis + NYnyr 2U+ n(n—1)Yn 
— (2n —1) yng. 2 —1(2n —1) Yn 
+ 1Yn 
BYnist Yn KM, 


CHAPTER XIII 


DIFFERENTIAL CALCULUS. EXPANSIONS. 
TAYLOR’S AND MACLAURIN’S THEOREMS 


1. It is often necessary to expand f(z) in a series of ascending 
powers of z This has been done by ordinary algebraic or trigono- 
metrical methods in such cases as (a + a)”, e”, log, (1 + x), sin a, etc. 

The various methods which can be employed may be summarised 
as follows: 

I. By algebraic or trigonometrical methods. 

II. By the use of Taylor’s or Maclaurin’s Theorem. 

III. By the use of a differential equation. 

IV. By differentiating or integrating a known series. [For 

integration, see Chapters X VI-XX.] 


2. I. An example will make this method sufficiently clear. 


Kaample 1. I ee = 
tt ea a 
1 ie 
71 oR wo 


by actual division. 


3. IL. Taylor's and Maclaurin’s Theorems. 
Assume that f(a +h) can be expanded in a series of positive 
integral powers of h. 
Let fa@th)=a+bh+ch?+dhi+.... 
Differentiating with regard to h, 
ff (@+h)=b4+2ch+3dh? +..., 
Sf’ (@th)= 2% +3.2dh+..., 
and so on. 
Put h=0, then we have ‘ 
PO? vases 


a=f(e), b=f'(@), o=15! 


Hence i ? . 
flath)y=fa)t hf (@)t sf’ @)+ gf" @) +. FIM @) 400. 


This result is known as Taylor's Theorem. 
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4. Ifin the above result we put «=0 we obtain 
F(h) =F (0) + hf" (0) + 5,F" (0) + bs +f (0)4 Re 
which, i. hel notation, may be written 
F(o)=f(0) + af (0) + FF" (0) + + ZF (0) 4. 


This result is known as Stirling’s or Maclaurin’s Theorem. 
The student should verify the common algebraic and trigono- 
metrical expansions by means of the above theorems. 


Example 2. . 
F (2) =log (1 + €), 3 f (0) = log 2, 
re ee ee +1) = 
¥ (@)= 747! 1 +e?’ f’ (0) =4, 
Ud e* id 
4 ©) - Gye f° O)=b 
m” oe ate , my -_ 
f° @)= ye f= 0, 
iv — ae re iv a S. 
etc. 
Whence, by using Maclaurin’s Theorem, 


of 


5. The above results are not universal in their application. Thus 
it can be shown that the method fails if f(«) or one of its derivatives 
becomes infinite or discontinuous between the specified range of 
values of z. Further, the series obtained must be a convergent series. 
Lagrange has shown that the remainder after the first n terms have 


n 
been taken from Taylor’s series can be expressed as vy f™ (a + Oh), 
nr: 
where @ is a positive proper fraction. The corresponding value of 
1) 
the remainder in the case of Maclaurin’s Theorem is = f™ (Oz). 


Unless therefore those expressions tend to vanish when n becomes 
infinite, the series is divergent and the method fails. 


EXPANSIONS 


6. III. The use of a differential equation. 
The following are examples of the use of this Beatie. 
Eaample 3. 
ae + 2)"=h+ 0+ 20+ a,0+.... 
Then =n(1+2)" or (1+ 2) y,=ny. 
Diderentatog the first equation, 


Yy = a, + 2a,” + 3a,0? +. 
Therefore 


(1 + &) (a, + 2age + 8a;,2? + ...) =n (AQ) + a2 + Aga? +...). 


Whence, by comparing coefficients, 
Oh = Nig, 
242 + 4, = NQ,, 
30, 4- 20, = Rasy; 
ete. 


Putting z=0 in the original equation we have 


a, = I, 
and successively os 
AQ, = NAy =n, 
n—1l n(n—1 
gata oP O=1) 
_ (%—2) | _n(n—1)(n—2) 
ag le ee 
etc. 
Hence (L+ay=1+ne+(S)at+(3) att... 
Example 4. y =sin # = a) + a," + a,0°+..., 
fein 
V1—a? 


and (1 —a*) y?=1. 
Differentiating (1 — a’) 2y%y,— 2ay,?= 0 


or (1 — a) y, = ay,. 
Now sin-! = — sin (— 2). 
Therefore 
Ay + Ay + A202 + a0 +... = — (My — Ge + a,2? — ase? + ...). 


Whence d=t,=a,=...=0. 
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Thus we may write 


Y=AT+ AgO+... + nyo + ..., 


Y — a, + Age? +... + Gong 2" +..., 


x 
Yi=Q, + 8a,0?+...+2n4+ 1 deny 0" 4+..., 
Y= 3.2asa+...+2n+1.2N. deny v"+.... 
Therefore 


(1—2")(3.2.aga+... + 2n 41.21. dongs 0" + ...)= 2(a, + Baza? +... 


+ 2n+ 1 engi 0" + :. F. 
Whence & = 3 . 2d, 


3a,=5.4a;—3. 24s, 


(2n + 1) dens, = (2n + 8) (20 + 2) Gongs — (20 + 1) 2n. Gongs 
or (2n +1) deny, = (2n + 8) (2n + 2) ants. 


But in the limit when « — 0, = 1. 


Therefore a, = 1, 
; 1 
and accordingly Ws= 5-5» 
1.3 
Se ab’ 
etc. 
; 12 1.8 2 
—l = = Ss ee .. 
and sin w= e+ 5.5 + aR te 


7. IV. Differentiation or integration of a known serves. 
The method is sufficiently indicated by the following examples. 
ge # 


Example 5. sInw=H— te, 
dsin a oF 
a =cos@=1—sotqi- 


We have therefore obtained the expansion of cosa in terms of « 
from that of sin a. 
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Example 6. Let 
y = log, (1 —#) =a, +a,% + a,27+..., 


dy il 
da Lag 1 t At + 80g? + Shee e 
But papaltatottaty.... (a <1) 


Therefore, comparing coefficients, 
a4=—-1, a,=—-}4, a=—}, ete. 


Also putting «=0 in the original equation, we find a,=0. 


Hence log. (1 —#) = —a@— gy — Zoe. (a< 1) 


in terms of # enables 


Thus our knowledge of the expansion of i Z = 
us to find the expansion of log, (1 — 2). 


H. 1.B,I. 6 


CHAPTER XIV 


DIFFERENTIAL CALCULUS. MISCELLANEOUS 
APPLICATIONS 


LIMITING VALUES OR UNDETERMINED FORMS 
1. In certain cases the value of a function or yee be at once 


ascertained. A common form is illustrated by 3 Li a = , the value 


of which is not apparent since it takes the form 7 The function 


in this case is said to take an Undetermined or Indeterminate Form, 
and since, in the limit, its value is represented by the ratio of two 
indefinitely small quantities, it may be expected that the methods 
of the Differential Calculus will enable the Limiting Value of the 
function to be found. 


2. The process to be employed can best be illustrated by examples. 


Example. To find Lt [== | 
Let 2 =1 +A, then 


u[2t]-» (Ga 
a baat 2) he + & aM | 
= 


h>0 


= Lt [n+(3 Vat (Rit... | 


=n. 


Eeample 2. To find Lt ee | 


=Lt E —t+50-. J (by actual division) 


since the series obtained by division is obviously convergent. 
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Example 3. To find if 9, 
20g” —I 


me ¢sj;,9t = 


. 4d [9” loge-+¥ (loge)?+... es 1] 


x2—>0 Ca 


2 
g Ke loge +5; (log ¢)* + Hog g 
2 2 l 2 
+ |eloge+F (loge) +...| — YS 
x 
wlogg +5 (log g) +... 


te slog glogc+termsinvolving zand higher powersof « 
i ne a. 5S a aeRO 
fee wlogg +5 (log g)' +... 

— Ti glog glogc+terms involving w and higher powers of « 


xr—>0 


log g + 5; (log g)? + es. 
=g loge. 
3. The formulas of the Differential Calculus can be used for the 
solution of problems of this class, in the following way. | 
Let J() take the form 2 when 2a. 


$ (#) 0 
Then if c=a+h, by Taylor’s Theorem, 


f(e) _ f(a) +f’ (a)+... 


$(a) ~ (a) +hg (a) +... 
Therefore, when «>a, and consequently h, f(a) and ¢(a) all 
—(, we have 
tp £@) _ £2). 
amap(t) (a) 


If this is still an undetermined form, the process can be repeated. 
Taking Example 3 above, we have 


f(@)=9" -9, (x) =g*—1, 
f'(aj=9" Ploggloge, (x) =g* logg. 


6—2 
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c# log g log c 
Lt PO ay Lo eagles « 
ro0P(t) 2a g*logg 
= g log c, as before. 


Therefore 


4. Other Undetermined Forms are found when a function in 
the limit assumes one of the following forms: 


00 
Oxo, —, w-—c, & wo of I”. 
oo 


These can all be made to depend upon the form >: 


Thus if f(x) and ¢(«) be the two functions involved, in the first _ 


ease by writing f(x) x (x) as [ te = we have the form ; and can 


proceed as before. 
Similarly for the second form. 


The third form can be written as 0 {1 — =I ; the second term 


within the bracket can then be evaluated, and the quantity within 
the bracket will then be found to be zero or otherwise. In the 
former event the product is of the form o x 0 and can be evaluated ; 
in the latter event the value of the product is clearly infinite. 

The last three cases are solved by taking logs, when the logarithm 
of the function will be found to fall within one of the cases already 


discussed. 


5. In connection with the form 1%, it is useful to remember 
that, from the development of the Exponential Theorem in algebra, 
it follows that 


Lt (1+ =) =¢ 
200 Yo 


and Lt (1 +4 =) = ¢*, 
x 


C3 00 


MAXIMA AND MINIMA 


6. If f(«) continuously increases with x until it attains a certain 
value a and subsequently decreases, the value a is said to be a 
maximum value of the function. Conversely if the value of the 
function continuously decreases, while «# increases, until it attains 
a value b and subsequently increases, the value 6 is said to be a 
minimum value. 


v 
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By this definition there may be several “maximum” values and 
several “minimum” values; further, a “minimum” value may con- 
ceivably be greater than one of the “maximum” values. 

These conceptions are illustrated in the following diagram. 

Thus the points A and D are maxima, whilst the points B and 

.& are minima. At the point C 
the value of the function ceases 
momentarily to increase but is in- 
creasing both before and after the 
point is reached. Such points, 

_ where the rate of increase or de- 
crease in the value of the function tends to become constant as the 
point is reached, are called points of inflemon. 


7. In obtaining criteria for ascertaining the maximum or 
minimum values the following considerations are of importance. 
If f(z) increases in value with 2, it 1s clear from the definition 


that a is positive. Conversely if f(x) decreases with 2, wy 18 
negative. 


Further, if f(«) is continuous and assumes two values f(a) and 
f (6), one of which is positive and the other negative, it follows that 
for some value of between the values a-and b, f(#) must be either 
zero or infinite. 


8. Now for a maximum value the function increases up to that 
value and then begins to decrease. Therefore at a maximum value 


dy must change from positive to negative; similarly at a minimum 


dx 


value 4 is changing from negative to positive. It follows from the 
preceding article that at a maximum or minimum point the 


value of = as zero or infinite. 


dy 


9. Further, at a maximum point 5 is decreasing from positive 
2 
to negative, and therefore = must be negative or zero. At a 


ug ay . oe 
minimum point the converse is the case and ed 1s positive or zero. 
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Hence in order to find the maximum and minimum values of 
J (@) it is necessary 

(1) to find the values of # which satisfy the equation f’ (x) =0 
or f’ (&) =ieo., 

(ii) to ascertain for each of these values whether f” (a) is negative 
or positive. . 

If f” (z) =0, a further test is necessary. For the investigation 
of these cases, which include values which give rise to points of 
inflexion, the student is referred to more advanced treatises on the 
subject. 


Example 1. What fraction exceeds its pth power by the greatest 
number possible ? 
Let « be the fraction; we have to find the maximum value of 


y= x2 — xP, 
dy _ - 
as = il — mF a 
For a maximum or minimum value oy = 0, therefore 
1 — pa? = 0, 
and gen = 1 
Pp 
or ma Pl, 
2 
Now oo —p(p—1)2?- 
_P=2 
=—P(p— lye? 


and is therefore negative. 
Hence the above value of x is a maximum value, 


Example 2. What are the dimensions of the largest rectangular 
box on a square base the area of whose surface does not exceed 
12 square feet ? 

Let x be the volume of the box, a the length of the base and 6 
the height. 

Then the surface is 2a? + 4ab = 12, whence 

6—a 


as 2a 
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Also x=arb 
— 772 
a? (6 — a’) ce ef. 
2a 
de _ 3 _ 30° 
da 2’ 
aa 
da? = — 3a. 
For a maximum or minimum value 2 = (0, therefore 
3a? 
3— Sos 0 
or a = 2, whence b = V2, 
Gx . 


and since —, is negative when a= /9, this value gives a maximum 


value for 2. 


CHAPTER XV 


RELATION OF DIFFERENTIAL CALCULUS TO 


FINITE DIFFERENCES 
1. By Taylor’s Theorem we have 


fletl=f@+f'@ +O... 


But f(«+1) = Ef(x), therefore if the operation = 


D we may write symbolically 


Ef(2)=(1+D+D?+...)f(2) 


be denoted by 


=e? f (a) 
or E=e?, 
whence A=eP-1 
and D = log, (1 + A). 
_ ar» A? & 
Therefore Cs A- 7 as — 12. | eee ee (1). 
If the interval of differencing be h, we have 
/ Re oy 
f(ath)=f(a)+hf' (a) + ait (a) +... 
and therefore A= ¢e*) —}, , 
Further, D? = [log, (1 + A)F, 
d?y ay A? = AS 2 
whence 7 [4-S+g7-| y 
= A?y — Aty + Atty —§Ary +... ...(2). 
Bindiigaly EY = Noy — BAY $ EDY Hoe seeeeeeeeee (3). 


2. The above formula (1) gives a convenient expression for the 
differential coefficient in terms of advancing differences. But from 
the nature of the differential coefficient, it may be expected that a 
better result will be obtained by applying the formulas of central 


differences. 
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Thus, taking Stirling’s formula, we have 
Af(O)+Af(—1). # 
f(a) = (0) +o SOF OIOD 4S arp) 


x (x®—1) A*f(—1)+ A f(— 2) 
a elle calle a 
Differentiating with regard to «, 


p(y O+APED _ AAC+ AC) , 


+ terms involving x and its powers. 
Putting «= 0 we obtain 


Sg oe ae .. (4): 
Taking the first term - we have 
f'O= Af (9) nie 1) -/E US 1) (5), 


a very useful approximation, which, by altering the origin and the 
unit of measurement, may be expressed more generally as 


ig) Leth =f(e-) 
f@- eae 


If h = } we arrive at the approximate result 


GOGO)... sccnerescorgoneersnn ct (7) 
=A (+A)? f@) 
= [A — $A? + BA — FAS + HA’ — ...] f(*). 
Comparing this with the value given by formula (1) we see 
that the error involved in taking f’(x) = Af(a-—}4) is 
[— dhe A® + pe At — AAP +... J f(a), 
which is approximately equal to — 4, A? f(a — 3). 
Hence we may write 


f' (@) = Af (2 —B) — Atl (@ HB) cee erences (8). 


3. The introduction of the second term in equation (4) gives 
the useful approximation 


Fy’ (0) -_ 8 {fQ) =e a {f(2) —f (= 2)} — (9). 


The other formulas of central differences, when treated as above, 
give these and other expressions for f’ (x). The above, however, 
are the most accurate and useful for general use. 
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4. Using the central difference formula, we shall obtain for the 
differential coefficients of higher orders, 


f” 0) =A*f(-1) -—- WAYS (- 2) ee (10) 
md 7" MLCD+OFE®) 
mA? f— EB) MPTOKs occ ve esc ensenes et (11). 


5. An example of the working of the above formulas is now given. 


Example. Find the first three differential coefficients of log a, 
when # = 2. having given the values of log, 1°80, log, 1°90, .... 


ae AS AS AS 

‘587787 | :054067 | — 000271 | — 000038 | 000607 

641854 | 051293 | — 000233 | — ‘000031 | ‘000005 
-| 693147 | 048790 : 000202 | — 000026 | 000004 

"741937 | 046520 : 000176 | — 000022 | :000004 

°788457 |. 044452 , 000154 | — 000018 

832909 | ‘042560 ‘ 000136 


875469 | 040822 
916291 | 039220 
955511 


Using the advancing difference formula (1), starting at the term 
2'00, and remembering that the interval of differencing is not unity 
but ‘1, we get 


f' (#) =10 (-048790 + 


002270 a. ‘000202 ‘000026 é. ) 


2 3 4 5 
= 50000. 
Similarly 
J” (e) = 100 (— 002270 — 000202 — 43 x 000026 — § x 000004) 
= — ‘2499, 
J” (2) = 1000 (000202 + $ x 000026 + § x 000004) 
= 248, 
The true values are, of course, given by the differential coefficients 
: ; 1 L 2 
of logx, which are respectively = gular When #=2 these 


become 4, —}, 4; the closeness of the above approximations is 
apparent. 

Formulas (6) and (9) give respectively for f’ (x) values of ‘50042 
and ‘50000. 


CHAPTER XVI 


INTEGRAL CALCULUS. DEFINITIONS 
AND ILLUSTRATIONS 


1. In Finite Differences we discussed under the heading of Finite 
Integration the problem of finding the value of f(0)+/(1)+... 
+f(n—1), or, changing the origin and the unit of measurement 
the more general series 


h[ f(a)+f(ath)+f(a+2h)+...4+f(a+n— Ih). 


The Integral Calculus deals with the value which this summa- 
tion assumes when / becomes indefinitely small. 


3 


2. Now let F(z) be a function such that f(z) is its differential 
coefficient. Then by definition 
fos F(a+h)— F(a) 
h>0 h 


or f(a)= —— +a, 


where aq, is a quantity that vanishes when h->0. 
Then we have 
hf (a) =F(at+h) —F(a) +ha, 
hf(at+h) =F(a+2h)—F(a+h) + ha, 
hf(a+2h)=F (a+3h)—-F(a+2h) +ha,, 


? 


hf(a+n—1h)=F(a+nh)—F(a+n—1h) + hap. 
By addition, 


AL f(@)+f(ath) +f(at2h)+...+f(a+n—Ih)] 
=P (a+ nh)— F(a) +h [a + a+... + an). 
Now if a denote the greatest of the quantities a,, @,... an, the 
last term is clearly less than nha, or (b — a) a, if a+ nh be put equal 


to b. This term therefore vanishes in the limit and we are left with 
the relation 


Leh[f(a) + fat h) +... +fO-h)] =F) - F(a). 
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This result is usually expressed in the form 


[ f@)do=FO)-F(o, 


where F’(«) is a quantity such that f(z) is its differential coefficient 
and dx represents, in the limit, the indefinitely small interval 
between the terms which are summed. 


b 
The expression i f(x) daz is called the definite integral of f(x) 


with regard to 2, and b and a are called respectively the superior 
and inferior limits of integration. 

Where no limits to the summation are expressed and we are 
concerned merely with the form of the function, we obtain the in- 
definite integral. 


3. A geometrical illustration of the process of integration may 
be obtained as follows. 
Let the equation y=f(z) be represented by the curve in the 
diagram. Let OA =a and OB=b, so that 
MA = f(a) and NB= f(b). N 
Further, let AB be divided into n parts 
each equal to h. Then the sum of the rect- fn 


=b 
angles shown 1s clearly equal to = h f («), and, 


in the limit, when h>0 the value of the 9 ee 
integral becomes equal to the area between the curve and the axis 


of x bounded by the ordinates MA and VB. 


CHAPTER XVII 


INTEGRAL CALCULUS. STANDARD FORMS 


1. In the preceding chapter we found that the process of integra- 
tion is the converse of that of differentiation. In other words, that 
given f’ (a) we have to find a function f(x) such that f’ (a) is its 
first differential coefficient. The analogy with the process of finite 
integration is apparent and the remarks made in Chapter VII, 
§5, apply equally to the present case. 

It may, therefore, be said again that the process of integration 
cannot be carried out for every function, that definite rules cannot 
be laid down to apply in every case, and that the student must look 
to applying the results of the differential calculus as a guide in 
solving any problem presented by the function under consideration. 


2. In the first place it is necessary to point out that the ordinary 
algebraic laws apply to the integrating symbol fd« in the same way 
as they have been shown to apply to the symbolic operations A and 
d 
aes . Thus: 

(1) The operation is distributive for, if u,v, w,... be any func- 
tions of a, 

rales + judx + furde + a sis 
and therefore 


[ude + nde + [wde+...= [(u+o+w+ jae ...CL). 


(2) The operation is commutative with regard to constants, for, 


p av d(cv) cdv_ 
ton 2a i 
therefore i cuda = cv=c | FS io... oos. eon (2). 
3. It should be added that the process of integration introduces 
| hae nel 
a constant into every indefinite integral, for, if w= ree 


d(vtc) dv _ 


a ie 
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and therefore I udx=VU+6, 


which result is, indeed, obvious from the consideration that a con- 
stant term disappears on differentiation. 


4. In accordance with the symbolic notation which has been 
previously developed (Chapter XV, § 1), we may write 


[7 («) dee = (s) i re. (8) 
= DF @) 
= Feces] ee (4). 
5. Standard Forms. 
: ce (n-+1)a" 
(a) @. Since SEL M402 _ on 


N+1 
therefore | a"dxz= z +c. 
n+ 


1 
. d(logz) 1 
1 EOL a 
(b) a. Since Ta = 
therefore | . =loga +c. 
—" d(e*) _ 
(c) e*. Since 7; e*, 


therefore | eda =e +. 


a” 
(@ &. Sine i srl) a ade 
igna a —t—(‘ié«i 
therefore | ada = a + 
log. a 


Other integrals may be obtained from a consideration of the table 
of standard forms of differential coefficients given in Chapter XI, 
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§ 3. The results are embodied in the following table. Their verifi- 
cation is left as an exercise to the student. 


Function Integral 


gut 
n+1 


[008 xdx 
fain xvdaz 


| sec? rda 


| cosec? 7 dx 


In each case the constant c has been omitted for the sake of 
simplicity, but its importance in the result must not be forgotten. 


6. Sums of the above functions can be integrated by the use of 
the distributive property of the operation of integration. Thus 
|e +- 2) de= |a% de + [2 de 
Wha 
log, 2 
Similarly | (ast + ba~*) dx =a | xdax +b i das 
“® ig 


= ee 


5 3 


+¢. 


7. The problem of integration may be presented in another way. 
Thus to find the sum of 
1 1 1 i 
a atm + e+2m "aba 
when « is indefinitely increased. 
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The series may be written as 


Lig 1 1 1 m 
oll wm > se Poe cal e 
re ae ag l+2.— 


The terms in the denominator consist of values of a quantity, which 
can be represented by y, increasing by equal increments of = . Also 


the sum of the series is multiplied by ~ , which is the increment 


in the value of the denominator. 
Since the initial value of the denominator y is unity and the 
final value is 1 +m, from our definition of an integral we may write, 


m 
when x->0o0 and therefore - 


1+m dq 
Lt a ewe here ie d jet = 


zo m | 1 142 142” L+e.™ zm y 


CHAPTER XVIII 
INTEGRAL CALCULUS. METHODS OF INTEGRATION 


I. THE METHOD OF SUBSTITUTION 

1. Various methods are available for treating the integration of 
functions which are not among the standard forms. Of these methods, 
one of the most important and most easily applied is the method of 
substitution. By this method the independent variable « is changed 
to y where y is a function of x. Thus an integral in terms of y can be 
obtained which can frequently be made to assume a standard form. 

For suppose y= F(a). 
Then dy = F’ (a) dx. 


So that if the original integral is of the form i SF [F (2) F' (2)dz, 


it can be expressed as { JF (y) dy, which may be a standard form. 


No general rules can be given, but the following examples will 
indicate methods which can frequently be employed with advantage. 


2. (1) |e te ba)” dx. 
Put a+ be= y. 
a bde = dy. 
at n+ 
Therefore [(atbeyrdem [Pay= yg tom Grae te 
aida 
} (a + ba)" 
Put Qa b wie Y. 
Therefore 
de =e (y— 2 dy 
(a (a+ bx)” 


Bo. Zay-" + q@? yr”) dy 
Se) aS hey FEY nl Se Be, 
«BL (n= 8)y™>— (n—-2) yy? (n-1I)y™ : 
where y= a + ba. 
M. T. 8.1 7 
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(3) [@+ 200) (a+b0+ ex)" der 
Then (b + 2cx) da = dy. 
Therefore 

Jo + 2cx) (a+ ba + caz*)" da = [yrdy => +k 
A (a + bax 4 cae)nh 
7” oh. 

(4) | { (log a dx . 
Put log a= y. 
Then Pain 
Therefore 

(log i de _[ ng (og aye» i 
| =[y y=syite= n+l 
d 

(5) fe 
Put J («)=y. 
Then J  («) da = dy. 
Therefore 

Cine WY = logy + o=log [f(a)] +0 
ada 

(6) T+? 
Put m= ¥. 
Then 2ada = dy. 

da _ ydy _ | i 

Therefore ent] i - + (1-35) 4 


=1[fav-[555] 


Diogecy 


=$[y—logl+y]+e 
= 4[a*?—log1+27]+<e. 
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sin «da 
7) lavbees: 
Put a+bcosrx=y. 
Then — bsin «dx = dy. 
Therefore 
snadx A1fdy 1 _ _log(a+bcosz) 
\ ese 7 pe b er 
2 sin = Cos = 
2 2 
d | sec 5 da 
2 tan 5 
Put tan 5=y. 
Then 4 sec? 5 de=dy. 
Therefore cs “ i dy _ log y +e 
sin & y 


= log tan 5 +. 


Corollary. Since 
cos # = sin @ + x), 


da Tr 
Lo 7 5) 
COs 2 eEree ¢ 75 ane 
II. RATIONALISATION BY TRIGONOMETRICAL TRANSFORMATION 
3. The above methods can be extended by the use of trigono- 
metrical functions and the relations which exist between them. 


The following are examples: 


dx 
° lane 
Let % =a Sin y. 
Then dx =acos ydy, 
Va? — a = acos y. 
7—2 


and 
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Therefore epee wae =-| dy=yt+c 
= sin —+¢, 
a 
dx 
2 ——<_—— os 
@) Va? — a? 
Let L=Asec y. 
Then dx = asec y tan ydy, 
and Va? — a? =a tan y. 
ceca 
+4) 
So = = [ SE =log tan (ae 9) +e 
1+ tan g 
= log 7¢ 
1 —tan 
wala: +siny 
cos ¥ 
a? 
= log —_—_ Te 
Pe 
i (* + Va — = 
a 
dx 
3 {= ° 
(3) ae 
Let 2=a tan y. 
Then dz = asec’ ydy, 
and Va? + a? =a sec y. 
Therefore 


dx  fdy _ i 
aia oh = log tan (T+ yy +e 
=log(secy+tany)+c [See Ex. (2)] 


2 
(= hte 


= log 
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From the above integral, two sther: _ fm portant: fortis’ An be 
obtained. 


da 
4 ees 
(4) Va + 2ba + ca2 


This can be written as 


__ te ened) 
exo Ve {(cw + b)' + ac —B} {(ca +b)? + ac— By} 


Disregarding for the moment the constant multiplier, this can 
clearly be transformed into one of the forms given in Examples 


(1)—(3) above. 


Corollary. Since 
dx = a ewer 
Freserees; 5 Ea Een) Xe 
wes yeaa =2 sin, /2=4 +¢. 


q (pe — qb) 
Gy {Kemer | iia 


Va + 2ba + ca? Va + Qba + cx? 
q (dat 2ba + ca?) (po? [ 
Va + 2ba + cx? ae 


The value of the first integral is 2 (a + 2ba + cut)? and the second 
integral can be dealt with as in Example (4). 


4. General remarks in regard to the method of substitution. 
When the precise substitution which will enable the integral to 
be solved is not apparent the device illustrated in the following 
examples is often of assistance. 


(1) (a+ a) dx 


GElice 
The appropriate substitution 1s not apparent, so assume 
Vei+e= y’, 


where r may have any value, to be determined subsequently. 
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Then! |2 0) 91 aah 
and a * 2r y— dy. 
The integral thus becomes 
(y*" +a—c) 2ry*—dy 
oe) led 


Clearly this integral can be evaluated if y” = y""—, ie. if r=1, 
when we get 


pe Grtag 4 affa+ ae pba) 


(a —b) il il 
—t sS ve —Vc—b arent 


a—b y—Nc—b 
ale 
ve — "By +e et 


= Bi) + 


where y=Vaz+<. 


(a —p)V (a+ 2ba + cx") 
Let us take the linear function as likely to lead to simpler results 


and, as before, put 
c—p=y" 
so that da=ry—dy. 


Then the transformed integral is 


f ry" dy =| rdy 
y' Va+2b(y+p)tc(y*+ py JyVa+2b(y" + p)t+c(y" +p)? 


Obviously if r =— 1, the denominator reduces to the form 
Va’ + Qb’y + c'y? 
and the expression can be immediately integrated. 
dx 
(3) = pee 
(a + cx)? 


This example is a little more difficult than the two previous 
ones. As before, let «= me: so that dz=ry" dy. Then we have 


ry" dy : 
(ei (a+ ca)? nee (a+ cy)t 
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The expression within the brackets in the denominator is of an 
even degree in y, so we will express the integral in the form 

| ry? d(y’) _ | ry’ dy?) 
2(atoy")! J 2y' (ay + 0)" 
This can clearly be integrated if r—2=38r, when r=—1. The 
integral then becomes 
-5| | e's d (ay? +c) _ 1 
2J(ayt+c)! 24) (ay®teo)t a(ay*t+e) 
a 
= ————__ +k. 
a(at ca)? 

5. Definite Integrals. The above examples of the application of 
the method of substitution all relate to indefinite integrals and the 
final form of the indefinite integral is expressed, in each case, in 
terms of x. 

In the case of a definite integral, a step may be saved in the 
operation by avoiding the final substitution of « in terms of y, 
provided appropriate changes are made in the values of the limits 
of integration. An example will make the point clear. 


2 wdxz 
Example. ae 
Put a = y. 
Then 2eda = dy. 
Also when « =0, y= 0 and when w=2, y = 4. 
Thus we get 


2 wde 4 dy Ta 4 ,/, 
rete Wwes-| +7 | = oat 
III. INTEGRATION BY PARTS 


6. From the Differential Calculus we have 
d(uv) dv, du 
dx da ast 

Whence, by integration, 


w= uw dx + [oS de, 
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dv du 
or [ug do=w - v dx. 


Consequently an integral of a function of the form u - can 


always be made to depend upon that of a function of the form v = 


The latter may prove to be a standard form and the use of the 
formula given above may enable an apparently intractable ex- 


pression to be integrated at once. 
The advantage of the method is best exhibited by applying it to 


a few elementary cases. 
(1) [ologeda= 5 | ogo 


=5 |e boge — [eB aa) 


d(2*) 9. 


a 1, ee 
(2) foe sf . da 


(3) | sin gdxz= {sin z da 


=osin2— fa A he 
i 


V1 
2 
=esin@ +5 st Sie =) 


7) er) 
=gxsin2+V1—a2+e. 
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(4) [Varta de = [Varro TF de 
=arva? +2? 


ere [x a aS 


wdax 
Va? + x 


2 
=Ve+e—[| Var | da 


ee a eee 2. 
ao arya [Varta de + [I 


C+ x2 


=aVa?+ a — 


Therefore 2 [Vat eda=aVata + a log [a + Va? + x] + ¢, 
and [Maree de =F Varta + Flog [a4 Varta + 5. 
(5) [log («+ Var a) da = flog (w+ Vat + at) $2 de 


=slog (04 Vee F a) — [eMC TNE TO) a 


=o log (e+ a4 + at) — - {| =. er 


=a log (2+ Va? + a?) — Va? + a? +. 


7. The above process can be expressed in a general form as 


follows: 
fucnat = | Ut —— dt 


= uw: fed - (F ({vat) dt. 


b 
To apply this formula to the definite integral | uvwdt, we can 
° a 


t 
write | v,dt as | v,dk, the upper limit being taken as ¢ since the 
a 


integral must itself be a function of ¢. 
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Thus we get 


b b t 
[ wedt= | mf maak — [Se (f ndk) at| 
b a b t 
=u | ade — ua [ mk f (I ndh) dt 


= [ ode [ S(f. mdk) dt 


since | ~ 0. 
Alternatively, if we express | vdt as the definite integral 


b 
_ | v,ak, we obtain 
t 


b ’) b b 
[ umdt=u | mdk + | Tit ([ ond) ae. 


8. Reduction Formulas. Where integration by parts is not im- 
mediately successful, continuation of the operation may ultimately 
lead to the evaluation of the integral. 


. a. da? 
Thus [rede = [ar dearer — fer de 


=v 2 [erode 
= giet —2 ¢ ate 
| dx 


= gie*—2 ve? — | e*da) 


= we” — Qre® + Qe* +. 


9. The above case is an example of a Reduction Formula, this 
designation being used since the application of the method of 
Integration by Parts effects the reduction by successive steps of 
one term in the integral. 


Thus 
- “ c gmti (1 uy a _ 
iE (l1—«#)"*dz= er jn a- x)" dx. 


* See J.I.A. Vol. 44, pp. 402-409. 


METHODS OF INTEGRATION 107 


Thus the degree of the term (1 —)" has been reduced by unity. 
Successive applications of the formula will reduce the degree 
further, until, if n is an integer, the term ultimately equals (1 — x) 
or unity. 

In the special case where the value of this integral is taken 
between the limits 0 and 1, we have 


1 m(J i ~ 1 gmt 1 1g 
fe (1 — 2) 22 |e - ay’ +f Sn —x)"" dx 


m+1 
nm fi 
=——— | a™ (1 — 2)" da, 
0 


whence by successive applications of the formula we find ulti- 
mately if n be an integer that 
nm! m! [amt da = n!m! 


m v7 — coe 
[v2 sora n+m! }o n+mtilt 


Many other important integrals can be dealt with in a similar 
way. 


IV. INTEGRATION BY THE USE OF PARTIAL FRACTIONS 


10. Any expression of the form where F(x) and f(a) are 


F() 
F(a)’ 
both rational integral algebraic functions of «, can be expressed as 
the sum of a number of terms of which the general forms are a,a” 
and amay & a For if the degree of #’(x) is equal to or greater than 
that of f(x), by division we can obtain a quotient, together with a 
new fraction in which the numerator is of lower degree than the 
denominator. The quotient provides the terms of the first form and 
the new fraction can be split up by partial fractions into terms of 
the second form. The integrals of these general forms are known 
and consequently the whole expression can be integrated. 

The following are examples: 

(1) [ PERS, - ! [cute _ ate de 


a /., b) a—bJ|a#—a «2—b 


“ea ~ [ag + p) log (w—a) — (bq + p) log (w@— b)] +. 


(2) | dx 2 dx 
av+2ba+e a Cease 
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The form of the integral depends upon whether ac —6? is positive 
or negative. 
In the former case the integral clearly is equal to 
a tan an + b 
Vac — b? Vac — b? 
In the latter case the integral becomes 


1 ’ az +b—Vb?—ac 
2 /b? —ac S as 5 be ee 


1 
(3) taro eae ,[e- "7 site| 
= | F — 80 — log (e+ 1) + 8log (w+ 2) | 


=8 log 3—9 log 2—§. 


CHAPTER XIX 


INTEGRAL CALCULUS. DEFINITE INTEGRALS. 
MISCELLANEOUS APPLICATIONS 


GENERAL PROPOSITIONS 


1. It is desirable to place on record several general propositions, 
in regard to change of limits, which are in the nature of being 
self-evident. 

In Chapter XVI it has been shown that 


[ s@de=FO-F, 


where f(x) is the differential coefficient of F' (x). 
It follows that 


b b 
1. i f(a) de = | POP cvcccccerccce (1), 
since neither x nor z occurs in the result. 
b a 
IL. | f(o)de=— i | ae (2). 
a 6 


Thus the interchange of the limits results in a change of sign of 
the definite integral. | 

On the left-hand side we have regarded the increment of dx as 
positive, so that, while # increases from a to b, the value of the 
integral is £'(b) — F(a). 

On the right-hand side the increment daw is negative and « 
decreases from b to a, giving a value for the integral of 


F(a)— F(b). 


UL. [ f@ae- [F@ dex + [s@de leer (3). 


For the left-hand side is F'(c)— F(a) and the right-hand side is 
F(c)—F (6) + £(6)— F(a). 


IV. [- f(a) de = | Sa ee (4). 
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For if we substitute a—z for z, then when «=a, z=0, and when 
o= 0, =a. 

Also —dz=dx. 

Therefore 


[ F@) dx =| -f@-2) d= ["f(a —z) a (by IT) 


. 


= | Je=nee Get 
0 
2. Differentiation of Definite Integrals. 


Let a | ) fide 


be a definite integral where the quantity c is independent of z, 
and the limits @ and b are independent of c. 


du 
To find a 


Ac inc. Then, since the limits are unaltered, 
b 
Au=| { f(a, c+ Ac) —f (a, c)} da. 


Therefore Au = | : Lf (@, ¢ + Ac) —f (a, ¢)} dx. 
Ac a Ac 


let Au be the change in wu corresponding to a change 


Proceeding to the limit, we have 
du [®df (a, c) 
a | ee ee (5). 
Thus the differential of the definite integral is reached by a 
process of differentiating under the sign of integration. 


3. An important use of this theorem is that of finding the values 
of other integrals from those of known form. 
For example, if the equation 


| e~ 8” dz = 
a, 


0 


* The student is referred to more advanced treatises for exceptions to this 
general result, 


MISCELLANEOUS APPLICATIONS Lo 


be differentiated » times with respect to a, we get 


n! 


iv) 
n pac el 


4. Areas of Curves. 

It has been shown that geometrically the definite integral 
| ‘" (x) da represents the area enclosed between the curve y= f(z) 
and the axis of 2 bounded by the ordinates «= a and # = b. 


Example 1. Prove that the area of the parabola 
y? = 4ac bounded by the curve, the axis of # and 


. ; ,ae—P 
any ordinate is two-thirds of the rectangle con- af 
tained by the ordinate and the intercept on the |, 
axis of x. 
Let ON = b. O 6b WN 


b a 
Then area oPN=| yds =| V4ax dx =2a* E oil’ 
0 0 
= 2(2a2 bt) 


=2PN.ON 
= # rectangle. 


5. Mean Value and Probability. 


Definite Integrals can be used to find the mean value of a 
function whose value is changing continuously by indefinitely 
small increments. 

Thus to find the mean value of f(#) for all values of x from a to b. 
If we divide b—a into n portions each equal to h, the mean value 


of the functions f(a), f(a +h),...f (a+n— Ih) is 
f(a)+f(ath)+... +f(at+n—Ih) 


n 


_A[f(a)+f(ath)+...+ f(at+n—Ih)] 
3 > b—a 

since nh = b—a. 

If now we make / indefinitely small, we shall have the mean 


value of all values of f(x) from f(a) to f(b). Consequently the 


— [r@ den 


mean value is 
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Example 2. A number nis divided at random into two parts; find 
the mean value of their product. 

Let x be one part, then the product is z (n—a). We have to 
find the mean value of this product for all values of x from 0 to n. 
Thus we get for the mean value 


1 {* 1/73 n? 
| e(n-2)de=> (FF) =F. 


6. The class of problem in Probability that admits of being 
treated by the methods of the Integral Calculus can best be 
illustrated by an example. 

Example 3. Three events A, B, C are known to have happened 
in the same century. What is the chance that the events happened 
in the order A, B, C; B happening within n years of the middle of 
the century? 

The chance that the events happened in the order A, B, C is 
clearly 4, if there is no limitation as to when B can happen, since 
there are six possible orders in which the events A, B, C can 
happen. 

To ascertain the further chance that, with the events occurring 
in that order, Bhappened within n years of the middle of the century, 
let « be the number of years from the beginning of the century 
to the event B,so that B occurs during the interval of time between 
aand«+dza. Then the A cases may have occurred in any of these 
x years, and the C' cases may have occurred in any of the following 
100 —a years. The total number of possible cases is therefore 
x (100 — x) in respect of the above interval of time da. 

But B must have occurred within n years of the middle of the 
century, therefore x can have any value between 50 —n and 50+7. 
Hence the required chance, being equal to the number of favour- 
able cases divided by the whole number, is 


| "(100 — a) de 
‘meaoo=nax (00) ~* (100) 


and the answer to the question is 


s[(ia)~*(05)] 30%) -8 0) 


MISCELLANEOUS APPLICATIONS 118 


7. An example of a somewhat different type is given in the 
following problem. 

Example 4. An event has always happened on an average once 
a year. Find the chance that it did not happen in a given year and 
prove it equals e7?. 

The event has happened on an average once a year. The chance 


i : . ' : 
that it occurred in any given nth part of a year is, therefore, —, and 
n 


the chance that it did not occur in that period is 1 ~*, 
By the ordinary rules of probability, the chance that it failed to 


occur in nm consecutive periods each of “th part of a year is 
n 
(1 --) . To obtain the solution to the problem, we proceed to the 


limit, thus obtaining Lt (1 -=)" , which by algebra equals e7. 


(ie aol) 
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CHAPTER XX 
APPROXIMATE INTEGRATION 


1. In many cases integration, or continuous summation of the 
values of a function, cannot be accomplished, either because the 
quantity to be integrated cannot be expressed as a mathematical 
function, or because the function itself is not capable of being 
integrated directly. 

In these cases formulas of approximation can be used, which may 
conveniently be divided into two classes, viz. formulas expressing 
the value of the definite integral in terms of 

(i) the sum of the successive values of the function and of its 
derivatives, 
or (11) the values of isolated values of the function, not necessarily 
successive. 

We shall now proceed to consider formulas of the first of these 
classes. 


2. The Euler-Maclaurin Expansion. 
Let Xf (a) = F(a), 
so that ” Sf (2) = AF (a). 
Then = f(@)=f(0)+f (1) +...+f(n—-1) 
=F(n)-—F(0). [See Chapter VII, § 4] 
Now F (x) = A f(a) 
=(e?—1)" f(a), since A=e?—1, 
which, by actual division, 
1 AnD 
- (5- 5 +g sagt 1.) f(o) 
oO = 
=[f@)de-46(@)+ kf @- ref" Ot 
Therefore 
F(x) F (0) =|" F@de-4{F@)-FO} +f FO) 
— rhe {Ff (nm) —F" (O)} +o - 
But F(n)—F(0)=/f(0)+ f(l)+...+f(m—1)). 
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Hence 


[F@ae=2fO+/Q)+-.+fO-Y +S @-H LS )-F'O) 
tbo (PM) —F"(O)}-.-. A). 

This result is not limited to the case where the ordinates are at 

unit distance apart, for, as has been remarked in connection with 


Finite Differences, by changing the origin and unit of measurement 
the formula can be given a more general form. Thus 


7 "f(@) d= tf(a)+fatr)t.. SNE) BG) 
ual we (a + nr) =e (a)} + 55 — ff" (a+ nr) —f" (a)} 


= Bai {Ff (at mr) — fl) te esesesereeceee (2). 


Example. Calculate 2 ; 


Taking values at unit Mitel we have, remembering that 


af. () 
ree and = mo 
dx #® = gl? 
Sepak 1, 1,1 1 1 
[oe 73 00tH tia tt Ba 35 Lv = + 30) 
1 
130 (~ a5 


+ 95) 


025000 
‘047619 
045455 
043478 
041667 
"020000 


223219 — 1, (0009) + +45 (00000369) 
223144, 


The value of the integral is clearly equal to log, 25, which to 
five places of decimals is ‘22314. 


3. Woolhouse’s Formula. 


This formula gives a relationship between the sum of consecutive 
ordinates and the sum of equidistant ordinates at greater intervals. 
It can be derived from the Euler-Maclaurin formula as follows. 


8—2 


116 INTEGRAL CALCULUS 
Formula (1) gives 
i Fe) da = Af Ot FQ) +A FO—V HAF) He FO) ~ 70) 
, + rho (F'n) —f'"(0)} - 
If the interval is = the formula becomes 
m{"sayde=4f(0+F (=) +F(=) + 
+£4f(0)- 3-1 @O- PO} +a I" -F'")} 
whence from the first of the above equations 
m |” f(e)de=m {f0) +f) +... +f} 
— 26 £0) +f) 3 SO) + ay LM) ~p"(0))- 
and from the second equation 
m |" f@)de= i +f (2) +4 FO} -2/O +L) 
— a Lf )-f'O)} + * Fe (f"(n) —f'"(0)} - 
whence, by subtraction, 
SO) = tet f(aam Sain» sats 
m= FO) +f} - 1 @)-7'O) 
+P @—f'"O}--. @) 


4. By putting a=0 and r=1 in formula (2) and proceeding as 
above, we obtain the following result: 


S(O) +f) +f(2)+. es n{f(0) si sae + 2 
+ f (mn)} — (if (0) +f (mn)} — ee (mn) — f’ (0)} 
we Cf" (mn)—f" (Pm veseee (4). 


This formula enables the sum of consecutive terms of a series 
to be expressed in terms of those at greater intervals. Thus the 
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work may often be shortened, particularly where the values of in- 
dividual terms of the series require to be obtained by calculation 
in the first instance. 


5. Lubbock’s Formula. 
Lubbock’s formula is similar to that of Woolhouse, but it sub- 
stitutes finite differences for the differential coefficients. 


Now D=A—$A?+4A8—fAt+ 
and D? = A’ — AS + 


Therefore, substituting in the above formula (3), we get 


FO) +f(=) +F(5) tt FO) 
= amide $f (n)} 2 FO) +f (n)} 


a 


(oy (n) — Af (0)} +5 — oy a (m) — A°F(0)} 


2 1)000=) {A%f(n) - a8 (0) 


sae SS eo Atf(O)}+... (5). 


6. By an alteration in the unit of measurement the formula may 
be expressed as 


AO AQ+ FO) Homan FO)+ FO) 2H) tS) 
mn)} "=" {Af (mn) Af (0)} 


+ wat {A*f (mn) — A*/(0)} 


— CADE =) tase (an) - af (0)} 


(n? — 1) (9n? —1) 

+ [asf (mn) — Asf(0)} +... 6) 

where A, A’, ... express the values of the differences taken over the 
interval n, and not those taken over unit intervals. 

In this form the formula corresponds to formula (4) given above. 
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7. In cases where the terms of a series decrease steadily and the 
term corresponding to the upper limit tends to the value zero, the 
above formulas can be simplified somewhat, since the final term of 
the series, its differences and derived functions all vanish. Also it 
is frequently the case that the end terms of the series are unim- 
portant and thus m and n may be taken at any suitable figures so 
as to correspond with tabulated values of the function. 

To illustrate this we will calculate, first by Woolhouse’s formula 


and then by Lubbock’s, the value of "ES (11)-2 
z=10 


Example 1. By Woolhouse’s formula. 

Let us use formula (4) and take n=15. The values of (1:1) 
tend to become unimportant when 2=100, therefore we will take 
m=6 and ignore the values of the differential coefficients at the 
upper limit. 


Then 
f(a) =(11), f’(#) =-(11)* log. 11, f’” (a) =— (1-1) (log. 11. 
Also. log. 1°1 = 09531. 
Thus we get 
& (t-1)-* 
10 38554 — (1°1)-10=2-6988 
25 09230 
2 os 
40 "02209 ae : x (11)! x 09531= °6859 
55 ‘00529 
70 00127 a x (11) “x (09531)°= 0234 
85 00030 
100 00007 


‘50686 x 15 = 7°6029 
And the final result is 
7-6029 — 26988 — “6859 + “0234 = 42416. 
The true value is, of course, 
n= 42410, 
a9 


which shows that the error involved by disregarding the final terms 
of the series is small. 
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Example 2. By Lubbock’s formula. 
In this example we will take n = 10, using formula (6). 


(1-1)-* A Aa 


38554 — ‘23690 +°14557 
"14864 09133 ‘05611 
05731 03522 02165 
02209 01357 

00852 

00328 

00127 

"00049 


one 99 x 1899 


62740 x 10=6° 720,000 
62740 x 10=6'2740 99. 899 


-_~. Bo 
480,000 “ one 


— 2°0238 


The result is 6°2740 — 20238 = 42502 as compared with the true 
value of 4:2410. 


8. The above examples illustrate points of disadvantage which 
may arise in connection with the formulas. 

Woolhouse’s formula can only be applied to a mathematical 
function where the differential coefficients can be obtained. In 
other cases, Lubbock’s formula must be used, and then it may often 
happen that the terms in the formula do not converge sufficiently 
rapidly to give a good result, unless a large number of terms is 
used. This accounts for the relatively poor result shown in the 
example above. If intervals of 15 had been used instead of 10, a 
worse result would have been shown. 

Some of these disadvantages can be met by the use of the 
formulas given in the following Articles. 


OTHER FORMULAS OF APPROXIMATE INTEGRATION 


9. The definite integral of any function can be expressed in terms 
of the individual values of any number of ordinates by assuming 
that the function can be represented, to a sufficient degree of 
approximation, by a parabolic function of the requisite degree in z. 


2 
10. Thus to express the value of | f (x) dx in terms of f(0), f(1) 
0 
and f (2), let f(a) =a + ba + cx’. 
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[N.B. Since we are to express the result in terms of three values 
of the function, we introduce three terms involving three unknown 
constants a, b and c in the expression representing f («).] 


2 2 
Then [ s(@)de=| (a+ bx + ca*) dx 
0 0 

ot Fadl 
= 2 Orie 
= 204 2b +e. 
Also f(0)=a, 
f)=at+bt+e, 


f(2)=a42b + 4e. 


These equations could be solved for a, 6 and c, and the resulting 
values substituted in the value of the integral found above. 


Alternatively, we may proceed as follows: 


2 
Let | f@de=pf)+afQ) +7F@). 
Then 20+ 2b + = patg(a+b-+ 0) +r(a+ 2 +40), 
so that pt+qtr=2, 
q+ 2r= 2, 
q + 4r = §. 
Whence r=4, q=4% p=4h, 
ant [" f@)de = LO2 YO FO) —— (7). 
20 


By analogy we may write 


| F@) dar = 5 {f(0)+4f (n)+f(2n)} vee, (8). 


This result is known as Simpson’s Rule. 


By dividing the whole area from 0 to 2n into consecutive 
equal spaces and applying Simpson’s rule to each, we obtain 


[ Fe) de= 51 FO)+ 4f 0) +2F@)+4fB)4 2 O+ 
+f (2n)}...(9). 
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This formula should normally yield better results than a single 
application of the formula (8) over the whole range of integration. 
It can, of course, only be applied where a sufficient number of 
values of the function to be integrated is available. 


11. Where the number of terms is four, by proceeding in a 
similar manner to the above, we arrive at the result 


[ ‘7 (e) do= 2 {f0)+8/(1) + 8F(2) + F@)} ---C10). 


12. The process of obtaining the desired formula may be sim- 
plified somewhat, when the given terms are arranged symmetrically 
about the central point, by adopting the central point as origin. 

This may be illustrated by proving the well-known Weddle’s 
Rule. 


6 
In this case | f (a) dx is to be expressed in terms of 
0 
FO), FL), «+» £6). 


Let f (a) =a ba + cx? + dat + eat + fa® + gx’. 
Then [-f@de= 6a+ 180+ — 


Also F (0) = a, 
FAY+ fC I) = 2a + 2c + 2 + 29, 
f (2) +f (— 2) = 2a + 8c + 826 + 1280, 
J (3) + f(— 8) = 2a + 18¢ + 162e + 14589. 
Solving these equations for a, c, e and g, and substituting the 


resulting values of these constants in the expression for the in- 
tegral, we obtain 


[°F @) de = z35 272f 0) + 17D 4FCD 
+ 216/(2) + f(— 2) + 417 (3) + f(— 3)}. 


This is not in a very convenient form for numerical work, so 
we add 
1 a, 
140 “Sf =—5t 207(0) + 15f(1) + f(- 1) 


— 6f(2) +f (—2) + f3) +F(- 3)}, 
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thus giving 
[7 F@) de + a5 MF (8) = Gay 252/04 427) FFE 
+2107@) + FCB) + 427 FB). 
Neglecting the term a A‘ f(—8) which will usually be very 
small, we arrive at the final result 


[7 @ do= Fa (f8)+5f2) +f) + 6f0) +f) 
- + 5f(— 2) +f(—3)}...CD), 
or, ill the origin, 
[. fieyde~ rp (f(0) + 5f(1) + £2) + 6F (8) +74) 
+ 5f(5)+/(6)}...(12). 
13. Another powerful formula can be obtained by expressing 


[t@ dx in terms of f (0), f(1), £(3), f(5) and f(6). 


Having five values of the function, let 
J (@) = @ + ba + cx? + da? + ex, 


Then [- f(a) da=6a+ 180+ “ore. 


Also Ff (0)=a, 
JF (2) + f (— 2) = 2a + 8c + 32e¢, 
Ff (8) +f (— 3) = 2a + 18¢ + 162¢. 


Solving for a, c and e, and substituting in the equation for the 
integral, we find 


| af (x) da = 2-2 f (0) + 1°62 { f(2) + f (— 2)} + 28 {f(8) + f(-3)}, 


or, altering the origin and limits, 
[£2 dr =n (2870) + Fn) +1627 @) + Fr) + 22/ Bn)} 


Similarly 


[ce de= n (287 Gay FRa) 41-625 Cay + FC) 
+ 2:2 f(9n)}. 
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Then | “F (a) de = | fa) d+ | —fO) Bie «x. 


= n [28 { f(0) + 2f(6n) +2. f(12n) +...} 
+162 { f(n) + f(5n) + f (Tn) + ...} 
+22 { f(8n) +f (9n) + ...}] 


Now in a series of values decreasing to zero, if 7n be so chosen 
as to fall just within or just without the limits of the table of the 
function to be integrated, we obtain the convenient formula 


} f («) da =n {28 f (0) + 162f(n) + 2-2 (8n) + 1:62f (5n) 
+°56 f (6n) + 1°62 f (7n)}...(14). 


The formulas in this Article are due to G. F. Hardy and the 
last is usually known among actuaries as 39 (a), since that is the 
number assigned to it in the original Text Book, Part II. 


14, A useful formula, involving only a simple summation of 
certain terms of the given series, is as follows: 


[Fe de= APL fea) +f (An) +f Gn) +f On). ..(15). 


By expanding f(#) in a series of ascending powers of « by 
Maclaurin’s Theorem, it is easy to show that the formula involves 
a small second difference error. It can, however, be used con- 
veniently where only a rough result is required. 


15. An alternative method of obtaining formulas of this char- 
- acter is to express f(a) in terms of the given values of the function. 
This can be done by Lagrange’s formula and the resulting ex- 
pression is then integrated between the desired limits. 

As an example, we will develop formula (18) in this way. 

The given values of the function are f(3), f(2), f(0), f(— 2) and 
f(—3). Then 


x u(2 x a+3)a(r2—2) (a2 
eos ee ee Chm as ce pai 
(a + 3) (w+ 2) (a@— 2) (#— 8) (a + 3) (x +2) a (a —3) 

Ett DE) 200) 4 Ore) 


4 (@ +3) (@+ 2) a (a — 2) 
~~ &xbxex1 NCP 
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Therefore 


[-f@de= f(-3) |" 
3 = ‘28 f (—3) + 1°62f(—2)+ 2-2f(0) + 1:62 f(2) + 28 f(3) 


i. (@ + 2) #(@—2)(#— 3) dx+ similar terms 
90 


as before. 


16. The above methods are perfectly general and afford means 
of expressing the value of an integral in terms of any given values 
of the function. 

It will be of value to indicate the method of application and to 
show the degree of accuracy in each case by applying them to an 
example. 


8 

Thus * _ 9 log, 2 = 1:38630. 
2 

(i) Simpson’s rule applied once gives 
$(4+4+4)= 142500. 


(i) Simpson’s rule applied three times over the values 2-4, 4-6, 
6-8 gives 

SS+EtE+ E+E +H + $) = 198770. 

(iii) The “three-eighths” rule gives 
GG +9494 4)= 140025, 

(iv) Weddle’s rule gives 

fs ($4+$4+444944454+4) =1:38679. 
(v) Formula (13) gives 

{28 +4) +162 (444) 4 5 = 138643. 


(vi) If formula (15) is used, n=°6 and we get 
lt. 2 ae ee ; 
4 Fare petra = 1:38839. 

As the values of the function are changing rapidly over the 
period used, the above is a somewhat severe test of the formulas. 
It will be noticed, however, that Weddle’s rule and formula (13) 
differ from the true value only in the fourth place of decimals. 


CHAPTER XxI 
PROBABILITY 


1. In considering the subject of Probability a clear appreciation 
must be obtained of the distinction between its mathematical or 
theoretical treatment and its arithmetical or practical development 
which is the basis of actuarial science. The difference between 
these two points of view will become clear ” examining a few 
simple illustrations. 

Probability or chance is merely an expression of relative degrees 
of uncertainty in relation to an event in respect of which our 
knowledge is incomplete. In the ordinary phenomena of life prob- 
abilities may not have a numerical value although distinctions 
are drawn. Thus we should say, in the absence of more precise 
knowledge, that a man of 30 was as likely as not to die before 
another man of 30, would most probably outlive a man of 60 and 
would almost certainly outlive a man of 80. The mind has clearly 
formed definite opinions on these points but the respective prob- 
abilities cannot be expressed in numerical form without further 
analysis. 


2. To obtain a measure of probability some unit must be taken, 
and for reasons of convenience it has been the universal custom to 
take absolute certainty as having unit probability. Starting from 
this standard, certain probabilities can be found from general 
reasoning, whereas others, such as the chance of dying in a year, 
or of becoming the father of twins, must be deduced from the ob- 
servation of suitable statistics. 


3. As an example of the former class, one may say that, on 
tossing a coin, it is certain that either “heads” or “tails” will appear, 
and as these possibilities are equally likely, the chance that heads 
will appear is clearly 4. 

We cannot however assume from this that, if we tose a coin four 
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times, two heads and two tails will necessarily appear. Indeed, as 
will be shown later, the respective probabilities are 
0 heads 4 tails 
1 head 83 tails 
2 heads 2 tails 
3 heads 1 tail 
4 heads 0 tails 


Jim [ae AR oko aR Oe 


What we understand is that, if the number of trials were 
sufficiently extended, the proportion of cases in which heads 
appeared would approximate more and more closely to $. It is in 
this sense that the mathematical definition of probability is to be 
understood in connection with actuarial work. 

The following Articles deal with the treatment of the subject 
from its mathematical aspect. 


4. From the foregoing considerations it will be evident that if an 
event can happen in ‘a’ ways and fail in ‘b’ ways, and each a these 
ways vs equally likely, the probability of ats happening ¢ 4s oon , and 
that of its failing is = - 5 For the sum of these two probabilities 


gives the chance of its either happening or failing, which is neces- 
sarily unity. 

Alternatively, we may say that the probability of an event is equal 
to the ratio of the number of cases favourable to the event, to the total 
number of cases. 

It follows that if p is the probability that an event = happen, 
the probability of its not happening is 1 —p. 


5. Another method of statement is to say that the odds are ‘a’ 
to ‘b’ in favour of the event or ‘b’ to ‘a’ against the event. 


Example. The chance of throwing a four at a single throw of a 
die is }, for there is one favourable result (a 4) and five unfavourable 
results (1, 2, 3, 5 or 6), all of which are equally likely. 

The chance of not throwing a four is §, or, in other words, the 
odds are 5 to 1 against the event. 
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6. The following proposition is also self-evident : 

If an event can happen in more than one way (all ways being, 
however, mutually exclusive), the probability of its happening at all 
is the sum of the several probabilities of its happening in the several 
ways. 

Thus if the chance of scoring a “bull” be +, that of scoring an 
“inner” be 3, that of scoring a “magpie” be 4, and that of 
scoring an “outer” be 4, the total chance of hitting the target at 
all (all events being mutually exclusive) must be 7, +44+4+4=2. 


7. The solution of elementary questions in probability depends 
therefore upon general reasoning, but calculation is aided in some 
cases by the theorems of permutations and combinations. 

Some elementary examples will now be given: 

(1) The odds in a given race against three horses are 11 to 4, 
13 to 3, and 7 to 2 respectively. Find the chance that one of them 
will win the race, a dead-heat being assumed to be impossible. 

The chance that the first horse should win is 54, 
, » second ” +; 
‘ » blamed : : 2. 
Thus the chance that one of them should win is 


ist ist $= 30- 

(2)..A has three shares in a lottery where there are three prizes 
and six blanks. B has one share in another, where there is one 
prize and two blanks. Show that A has a better chance of winning 
a prize than B in the ratio of 16 to 7. 

To be successful A may draw either 3, 2 or 1 prizes. 

He may draw 38 prizes in 1 way. 


He may draw 2 prizes and 1 blank in (3) x 9) = 18 ways. 


He may draw 1 prize and 2 blanks in et) x (3) = 45 ways. 


The total number of ways in which he can win at least one prize 
is therefore 1 + 18+ 45 = 64. 

Now three tickets can be selected in (3) = 84 ways. 

Therefore A’s chance of success is $= 4$. 

B’s chance is clearly 4. 

Therefore A’s chance is to B’s chance in the ratio 16:7. 


Alternatively A may draw all blanks in (3) =20 ways. His 
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chance of non-success is therefore 29 and his chance of success 
1—2¢= 48 as before. 


(8) If four cards be drawn from a pack, what is the chance that 
there will be one from each suit ? 
Four cards can be selected from the pack in i< ‘y= 270725 ways. 


Four cards can be selected so as to be one from each suit in 
13‘ = 28561 ways. 
Therefore the required chance is #8581. = 


(4) Out of a bag containing 12 balls, 5 are drawn and replaced, 
and afterwards 6 are drawn. Find the chance that exactly 3 balls 
were common to the two drawings. 

The total number of ways of making the second drawing is 


(G)= 924, 


To comply with the conditions, it must contain 3 balls out of 
the first 5 chosen, and 3 balls out of the 7 left on the first choice. 


The respective ways of making these selections are (3) and (A 
The total number of selections favourable to the event’is therefore 


(3) x (3) = 350, and the required chance is $59 = 28. 

(5) Twelve persons take their places at a round table. What is 
the chance of two particular persons sitting together ? 

Let the two persons be A and B. Then, since we are only con- 
cerned with the relative positions of the persons, we may regard 
A’s place as fixed. There are then 11 other seats in all, 2 of which 
are adjacent to A. B’s chance of occupying one of these is there- 


fore +. 


(6) What is the chance of throwing more than 10 in a single 
throw with two dice ? 
A score of more than 10 can be made by the following throws: 
5 and 6, 6and 5, 6 and 6. 


The total possible number of combinations is 6 x6=36 and, 
as 3 of these are favourable to the event, the required chance is 
go = 7s 

It should be noted that exactly 11 can be thrown in two ways, 
since (regarding the throws as consecutive) either a 6 or a 5 may 
appear first, but 12 can only be thrown in one way. 
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8. If the chance of an event is p and the measure of the 
quantity dependent on the event be a, the product pz is called the 
expectation. 

Thus we may have the expectation of a person in a lottery, or the 
expected value of a prize, or the expected number of deaths among 
a given number of persons, and so on. 

Examples. 

(1) What is the expectation of a person who is to draw one 
envelope from a bag which contains one £1 note, two 10s. notes and 
three blank pieces of paper, each placed in an envelope of uniform 
size? 

The chance of drawing the £1 note is 4. 

The chance of drawing a 10s. note is 4. 

The value of the expectation is therefore 

4 x 20s. + 4 x 10s. = 6s. 8d. ° 


(2) If the chance of dying in a year is s,, the — number 
_ of deaths among 100 people is 3 x 100=5. 


9. In the foregoing articles we have considered what are, in 
principle, single events. We have now to determine the appropriate 
formulas for combinations of two or more events. 

In considering problems of this kind, a close watch must be made 
to see if the events are dependent or independent. 

Thus, a bag contains 6 white and 4 black balls, and it is desired 
to estimate the combined chance of drawing at the first draw 
3 white balls and at the second draw 3 black balls. If the balls are 
replaced after the first draw, the second event is clearly independent 
of the first. But, if they are not replaced, the drawing of 3 white 
balls at the first draw will obviously affect the chance of drawing 
3 black balls at the second draw, and the two events will be 
dependent. 

Dependent events are also called contingent. 


10. We are thus led to the following proposition: 

The chance of two independent events happening is the product of 
the chances of their happening severally. 

For if the first event can happen in @ ways and fail in 6 ways, 
and the corresponding figures for the second event are a’ and 0’, the 
total possible combinations of events are (a + b)(a'+ 0’). Of these 


HT. BL ) 
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Both events may happen in aa’ ways. 
The first event may happen and the second event 
fail in ab’ ways. 
The first event may fail and the second event 
happen in ba’ ways. 
Both events may fail in bb’ ways. 
The al of their both happening is therefore ESN CET 
or = Sa , Le. the product of the respective chances of their 
happening. 
If the respective chances are p and q, 
the chance that both happen is PY 
the chance that the first happens and the 
second fails is p(l—4q), 
the chance that the first fails and the 
second happens is (l-p)q, 
the chance that both fail is (1 — p) (1 — q), 
the chance that at least one happensis 1—(1—yp)(1—q) 
the chance that one and only one happens 
is p(l—q)+(—p)q 
=pt+q—2pq. 


11. The above results can be applied by a slight modification 
of reasoning to events which are not independent. 

Thus if p be the chance of the event happening, and q be the 
chance of a second event happening when the first has happened, 
then the chance that both should happen is pq. 


_ 12. By successive applications of the above reasoning we can 
arrive at a formula for any number of events. 

Thus if p be the probability of an event A; and, when A has 
happened, q be the probability of another event B; and when A and 
B have happened, r be the probability of another event C; and so on 
for any number of events; the chance that all the events will happen 
as the product pqr.... 

Tt follows that if p be the chance that an event will happen in 
one trial, the chance that it will happen in each of a series of r 
trials is p”, and the chance that it will happen at least once in a 
series of r trials is 1 —(1—p)’. 
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13. This result can be developed more generally in the following 
way: 

If p be the probability of an event happening in one trial, what is 
the probability of its happening once, tunce, three times, ... exactly 
an n trials? 

Let q be the probability that the event does not happen, so that 
p+q=l1. Then if the event is to happen exactly r times, it must 
happen in each of a given combination of r of the trials and fail in 
each of the remaining (n—1) trials. The chance of this occurring 
is, as seen above, p’g"-". But the particular set of r trials can be 


chosen in (*) ways, each of which is equally likely. Therefore the 
total chance of its happening exactly r times is fe ) pq”, or the 


term containing p” in the expansion of (p + q)”. 

Thus the successive terms of this expansion represent the 
probabilities of the event occurring respectively n, n—1, n—2, ... 
times in n trials. 


14. The foregoing propositions can best be grasped by applying 
them to certain examples. 

(i) If four cards be drawn from a pack, what is the chance that 
there will be one from each suit? 

Let one card be drawn, which may be “ any suit. The — 
that a second card is of a different suit is 29, for there are 51 cards 
remaining, 39 of which will be of different suits from the first. 
Similarly the pn of drawing a third card of a different suit from 
the first two is 28, and that of a fourth different card is 43. 

The cagbieagal ey is therefore $9 x 26 x 48 = BASS. 

This result should be compared with that obtained in Example 3 
of §7, where the problem was treated as that of a simple prob- 
ability. It is clear that the same result must be obtained whether 
the cards are treated as being drawn simultaneously or successively. 

(ii) A man throws a six-faced die until he gets an ace; he is to 
receive £1 if he succeeds at the first throw, £4 if he succeeds at 
the second throw, £4 if he succeeds at the third throw, and so on; 
given that log. 6 = 1:79176, find the value of his expectation. 

The chance of succeeding at the first throw is 4. The chance of 
succeeding at the second throw is compounded of the chance of 
failing at the first throw and succeeding at the second throw, ete. 

9—2 
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The value of the expectation is therefore 
4+$8.4.$4+(8%.4.44... 
=4£{6+().3+(@*.3+...} 
on 5 i log. (1 - 8) 
= } log, 6 
= ‘3583 = 7s. 2d. 


(iii) Let it be assumed that the probabilities of dying within 
ten years after the ages specified are, on the average, as follows: 


Probability of dying 
within next 10 years 


What is the chance 

(a) that a person A now aged 30 should die between the 
ages of 50 and 60; 

(b) that two persons A and JB, aged respectively 30 and 
40, should be alive 10 years hence ; 

(c) that of two persons A and B, aged respectively 80 and 
40, A should die between the ages of 40 and 50 and B should 
survive to the age of 60? 

(a) The required chance is compounded of the chances that 
A should survive successively to ages 40 and 50 and should then 
die within the next ten years. 

The several chances are (1 —);), (1 — 4) and 4; and the required 
chance is the product of these three factors, namely 

HX § x3 = ies. 

(b) The respective chances of surviving ten years are (1 — 7) 
and (1—4). The chance that both A and B should survive that 
period is therefore 

Th X 9 = $7. 


(c) The chance that A should die between the ages of 40 and 50 
is 44 x 4. 


The chance that B should survive to the age of 60 1s § x 4. 
The required chance that these two events should both happen 


is, therefore, 
(14 x 4) x (§ x $) = ahis- 
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(iv) The faces of a die are marked with the consecutive numbers 
1, 2,...6. What is the chance that, after seven throws, the sum of 
the numbers exhibited equals 30 exactly ? 

This example is important, as it illustrates a method which can 
frequently be employed. 

The number of ways in which the seven numbers exhibited can 
total 30 is given by the coefficient of 2 in the expansion of 

(a+ 2+... + 2°), 
for this coefficient arises from the combination of the indices of a, 
taken together in such a way as to produce a, total of 30. 

Writing (v+a?+...+2°) as 2’ (1—a*)'(1—«)~ it is seen that 

we require the coefficient of x* in the expansion of 
(1 —af)' (1 —2)7 = (1 — Ta + 21a? — 35084...) (1-2). 

The expression within the first bracket need not be expanded 
further, since no term higher than 2” is required. It remains to 
combine the given terms with appropriate terms taken from the 
expansion of the second bracket in such a way that the power of 
x given by the product of the two terms is 23. 


Thus we get 
(1) x (Fae) = 475,020 o* 
Lx... X23 4.\ 
De ion, KG 
(214") x ( wos o)= 259,896 2% 
om.... K 11 
= 18 ‘a> ee OS ae = 23 
( 85a") x (Fo at) 16,170 a 
12,117 a 


The number of ways in which a total of exactly 30 can be made 
is thus 12,117. 

Now, on any of the seven throws, any of the six numbers may 
be exhibited. The total number of possible combinations of num- 
bers is, therefore, 6’, The number of combinations giving a total 


of 30 being 12,117, the required chance is a 
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15. The following miscellaneous examples are taken from the 
examination papers of the Institute ; they illustrate various devices 
which may be employed with advantage in the solution of questions 
of this character. 

(i) If x whole numbers be multiplied together, find the chance 
that the last digit of the product is a five. [1910, Paper J, Q. 7.] 

The chance is compounded of the separate chances that none of 
the final digits of the n numbers is even and that one at least isa five. 

The required chance is therefore 


‘ a 5” —_ 4” 
ar -@} =a 
(11) Four coins are tossed together and A is to receive £2 if 
exactly 2 heads turn up, and to pay £1 in any other event. Find 
Piers . the probability that after four trials A is £1 out of pocket. 
[1913, Paper I, Q. 8.] 
To satisfy the conditions A must win once and lose thrice. 
The chance that exactly 2 heads turn up equals the middle term 


in the expansion of (1 +4)'= a) x(4i= 2. 


Therefore the chance that he should win once and lose thrice 
equals the second term in the expansion of (§ + $)', which is 


4.3. (8) = foe 

(111) If a coin be tossed 15 times, what is the probability of 
getting heads exactly as many times in the first 10 throws as in 
the last 5? [1915, Paper I, Q. 9.] 

If n coins be tossed the chance that exactly r heads turn up 
equals the (r+ 1)th term in the expansion of (4 + 4)”. 

Now in the last 5 throws we may get 0, 1,... 5 heads and we have 
to combine the chance of any of these with the chance of getting 
the same number of heads in the first 10 throws. Thus we have 
for the several chances of getting 


0 heads in both sets of throws= (4). (4) 


lhead  , » = 54). 10)" 
2heads_,, -, » =10($)'. 45(4)" 
ZS oo» ml) OGD. Way 
4 ,, i, s 2 = 5(4)°.210(4)" 


5 ” ”? ” 9 = (4)° o 252 (4) 
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Therefore the total chance that the same number will turn up in 
both sets of throws is the sum of each of these distinct probabilities, 
3003 

(iv) A and B throw for a certain stake, each throwing with one 
die; A’s die is marked 2, 3, 4, 5, 6, 7, and B’s 1, 2, 3, 4, 5, 6, and 
equal throws divide the stake; prove that A’s expectation is 47 of 
the stake. : 

What will A’s expectation be if equal throws go for nothing? 
[1911, Paper I, Q. 8.] 

Let e equal A’s expectation. 

The chance of A’s throwing any of the numbers marked on the 
die is the same for each number. 

If A throws a 2, B must throw a 1 if A is to win the stake, or 
a 2 if A is to divide the stake. The chance that A throws a 2 is 4; 
the chance that he then wins the stake is } and that he divides the 
stake is also 4. His expectation if he throws a 2 is therefore §+4.4. 

If equal throws. go for nothing, A’s expectation after an equal 
throw clearly remains at e. In that case his expectation if he throws 
a2ist+de. 

We therefore have the following scheme: 


1.€. 


WY “ ' Chance | 4A’s expectation if equal | A’s expectation if equal 
S$ throw | thereof | throws divide the stake throws go for nothing 


$+3.3 : $t+te 
443-3 4+he 
$+3-4 st+he 
$+4-3 the 
ore. a“ 


The total expectation being the sum of the separate expectations, 
we have in the first case e = 4%, being the sum of the figures in the 
third column of the above table, and in the second case e = 7 + se. 

Whence e = 34. 

(v) A man tosses 20 pennies and removes all that fall head up; 
he then tosses the remainder and then removes all that fall head 
up, and so on. How many times ought he to be allowed to repeat 


this operation if he is to have an even chance of removing all the 
pennies before he has finished? [1906, Paper I, Q. 11.] 
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The problem is clearly the same if all the pennies are tossed 
each time; we then have to find the chance that all the pennies 
have turned up heads once at least. 

If n be the required number of throws, the chance that any 
particular penny has turned up heads once at least is {1 —(4)"}. 

The chance that all the pennies have turned up heads once at 
least is therefore {1 — (4)"}” and by the terms of the question this 
must equal 4. 


1 
Solving this equation we find (4)”=1-—- (4)°0, 
mn = 4°87. 


anny 


(vi) 2" players of equal skill enter for a tournament; they are 
drawn in pairs, and the winners of each round are drawn again for 
the next. Find the probability that two given competitors will play 
against each other in the course of the tournament. If n=5, show 
that the probability that a given player will either win or be 
beaten by the actual winner is ;. [1912, Paper I, Q. 6.] 

As regards the first part of the question, the total number of 
games that will be played is ) 


Qm14 M84 1 =I 1, 
Qn (21) 
a. We 


Also two players can be selected in ways. The chance 


that two will meet is therefore the quotient of these two values, 
Le. oni" 

The second part of the question can be proved by an inductive 
process. 

Let wu, be the required chance. 

Then in the first round he may either win or lose. If he wins, 
he passes into the next round and his chance of winning or being 


beaten by the ultimate winner then becomes w,p_,. If he loses, his 
opponent passes into the next round where the latter’s chance of 


, ‘ ‘ ws 1 
becoming the ultimate winner is increased to 5— 


Qn-1 ‘ 
We thus have 


Un = 5 Un + 5 + oni: 


PROBABILITY 4 137 


— 1| 1 mh * il 
Similarly 9 Un-1 = 7 Una a 4° yn=2? 
1 1 eh 
Fag er 
etc. etc. 
1 1 le 1 


ee a a 


As u, obviously equals 1, we have by addition 
1 n-1 n+l 


| ae oe 


It is instructive also to prove the first part of the question 
inductively. 

Let wu, be the required chance. This chance is compounded of 
the two chances that they meet in the first round, or that, not 
having met in the first round, they both pass into the second 
round where the corresponding chance becomes wy_). 
ay 1? the chance 


Therefore we have 


The chance that they meet in the first round is 


n 


that they do not meet is 


2 
* aiflgn gal Sink 
Po =e gee yy’ 
LS : 
or Ge ots 5 — 
n—1 __ ® n—2 __ 
Similarly (- 5 *) a =5 + Z Un—2 
Qn-2 _ ] 1 Qe] 
( 4 ) Un—o = 4 ar 8 Un—s> 
etc. etc 
23 — 1 il 2?—] 
( Qn—s ) a Qn—8 Qn—2 Us, 
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By addition 
ft @ i 1 
(27 — 1) tn=14+5 tt aa s+ TF onma + Sami 
(since u, obviously equals 1) 
1 
198 ony 
™ ae - ower 
19 
Whence Pe 
Un = Qn-1? 


as before. 


16. The application of the Integral Calculus to problems of 
mean value and probability is shown in Chapter XIX, § 5-7. 


EXAMPLES 
CHAPTER III 


1, Given 
f(1)=(@-2)(e-3), £@)=(@-7)(e-5), f(8)=(@-10) (e+), f(4)=68, 
obtain a value for 2, assuming second differences are constant. 
2. Find the zth term of the series 2, 12, 36, 98, 270, 768, etc. and the sum 
of m terms of it. 
3. Given that 
f(0)=66566, £(1)=65152, f(2)=63677, f(3)=62136 
find f(9). i by ae ? Ft (3) ? 
‘ 62% +29 
oer Ca (Be+4) (3x47) Ba+ sae 
5. Find the value of 
A” (32-+1) (82+4) (82+7) ... eat sa 
6. Prove that 
(n= 1)" f (2) 4m (n= Iy"“1 fle 1) + MD em —aym=2 fw 42) 
+...4+f(2+m)=n™ f (x)+mn™— Af (x) 
m(m—1 
RO) wt aif (e+... +f (@) 
Hence find the sum of the series 


+m. nee = 32+ ...+(m+1)% 


-: CHAPTER IV 
1. Given 
f(0)=70795, f(1)=72444, /(2)=74131, and £(6)=81283, 
find f(3), f(4), and f (5). 
2. Find f (35) given 
f(20)="01313, f(30)=-01727, f(40)=-02392, f(50)= 03493. 


3. Supply the missing term in the following table: 
J (0) =72795, 
f(1)=71651, 

f (2)=70458, 


J (4)=67919, 
T (5) = 66566, 
SF (6)=65152. 
4. Given 
f=, f(3)=18, f(6)=74, and f(9)=522, 
find the intermediate terms. 
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.5. Given 
f(0)=982038, (1)=97843, f(2)=97459, /(3)=97034, 

find 7 (2°25). 

6. Given 

Ff (0)=98023, ((10)=97651, /(20)=97246, £(30)=96802, 

find f (15). 

7. Given a: 
f (0)=58'842, f(2)=55°257, f(4)=51°368, (10) =37-977, 
complete the series f (0), f(1), ... £(10). 

8. If you were asked at very short notice to obtain approximate values for 


the complete series f(0), f(1), (2), ... £(20), being given that /£(0)=-013, 
J (10)='248, f(15)=°578, and £(20)=-983, what methods would you adopt, 


and what value would you obtain for f(9)? 
9. FQ)=15 FA)+F(B)=541; F(A) +f (5)+f(6)= 18-47 ; 
F(T) +F (8) +F (9) +F (10) +f (11) +f (12)=90°36. 
Find the value of f(x) for all values of x from 1 to 12 inclusive. 
10. Apply Lagrange’s formula to find f(5) and /(6), given that 
f(1)=2, f(2)=4, £(3)=8, f(4)—16 and f(7)=198; 
and explain why the results differ from those obtained by completing the series 
of powers of 2. 
11. Find the simplest algebraic expression in + which has the values 5, 3, 
9, 47 and 165 when z has the values 0, 1, 2, 3 and 4 respectively. 
12. Prove the following formulas for approximate interpolation: 
F(1)=f(8)- 3 [F(5)—F(-3)] 4-2 1F (-3)-F(- 5] eeeeeeees (1), 
FO) =4L/ 01) 4 f(-D)]- 3B (48)-F (}-3 FE1)-f(-3)}] «..(2), 
and apply them to find the logs of 45, 46, 47, 48, 49, being given 
log 42=1°62325, log 50=1-69897, 
log 44=1°64345, log 52=1-°71600. 


e 


13. Given 
10 10 10 
2f(#)=500426, 3f(r)=329240, 2 f(«)=175212 and f(10)=40365, 
4 
find f (1). 


CHAPTER V 


1. Use Gauss’ interpolation formula to obtain the value of f (41) given 
Ff (30)=3678°2, f(35)=2995'1, f(40)=2400'1, 
f(45)=1876-2, £(50)=1416°3. 
Verify your result by using Lagrange’s formula over the same figures. 
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2. Given the following table find (28) using Stirling’s formula: 
f(20)=98450, f(25)=96632, (30) =94472, 
f (35)=91852, f(40)=88613. 

3. Prove that 

2 2F (_ a 
F@)=3lf-H+sQraase pete. LCOtesed, 
and apply the formula to find f (32) given 
Jf (25)='2707, £(30)=°3027, /(35)=°3386, (40) =°3794. 


4. From the table of annual net premiums given below find the annual net 
premium at age 25 by means of Bessel’s formula: 


Age Annual Net Premiums 
20 01427 
24 01581 
28 ‘01772 
32 01996 


5. Use Everett’s interpolation formula to complete the series (25) to f(35), 


given that 
f(15)=305, f(20)=457, f(25)=568, £(30)=671, 


f(35)=897, f(40)=1190, f(45)=1481. 


CHAPTER VI 


1. Given the following table of f(z): 
f (0)=217, f(1)=140, Sf (2)=28, f(3)= —&, 
show how to find approximately the value of x for which the function is zero. 


2. Given that, when 
x=0, je (z)=0, 
z=1, f(2)=100, 
z=2, f(x)=2000, 
find # when f(%)=1900 by Lagrange’s formula of interpolation (applied 
inversely) and explain why the result does not agree with that found by using 
the formula f (7)=(1+ A)* f(0) and solving the quadratic. 


3. The following values of f (7) are given: 
f(10)=1754, £(15)=2648, /(20)=3564. 
Find, correct to one decimal place, the value of « for which f(x) =3000. 


4, f (30)=—30, £(34)=-13, £(38)=3, (42)=18. 
Apply Lagrange’s formula of interpolation inversely to find x, where f(x)=0. 
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CHAPTER VII 


l. Find the value of 
1 


and use the result to find the value of 
7 1 
2a (nel)(n+2)" 
2. Prove that 
3 (+1) olan. (n+). 


3. Show that 
% 2 
0% (2)= 25 | oe) 7 Ab (2) + ea Ath )-... | 


where ¢ () is any rational integral algebraic function of x. . 


4, Prove that if a, diminishes as z increases and converges to the limit zero, 
the sum to infinity of the series a,—a,+a3... is the same as the sum to infinity 
of the series 4a, —}Aa_+}A7a; — etc. 

Find the sum to infinity of the series ~y— 74+ 74... true to four decimal 
places, 


5. Prove that if the fourth and higher differences are ignored the sum of n 
successive terms of a function of which f(0) is the central term is 


nf (0) +5" atf(—1), 


where 2 is an odd number.. 


CHAPTER VIII 


1. Given the following values of f(z): 
sa f (1:41) =°7092, f(1°49)= 6711, 
J (1°52) =°6579, £(1°53)=°6536, 
find f (1°45). 


2. Use Divided Differences to find / (80) to the nearest integer, given 
f(70)=235, £(71)=256, £(79)=436, f(81)=484. 


3. Given 
f(20)='342, £(23)=391, £(31)='515, £(34)=-559, 
find f(30) by means of Divided Differences, and check the result by applying 
Lagrange’s interpolation formula. 
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CHAPTER IX 


1. The following values of f(a, y) are given: 
f (35, 55)=10°020, f(35, 50)=11:196, (35, 45)=12-019, 
f (40, 55)= 9°796, (40, 50)=10'894, (40, 45)=11-641, 
f (45, 55)= 9°583, (45, 50)=10591, (45, 45)=11°243. 

(i) Using only six of the above values, find f (42, 52). 

(ii) Making use of all the data calculate f(44, 51). 


2. Prove that if (0, 1)=/(1, 0) and f(0, 2)=/(2, 0), then 


Fe N=F0,0)+e49)] ant =S— ats |70, Obey 40, 1-7, 0) 
Find 7 (39, 33), given . 

f(35, 35)=3'151, (35, 45) =3°912, 

Ff (35, 40)=3-471, £(40, 40)=3-766, 


CHAPTER XI 


1 
1, Find the differential coefficients with regard to x of log xz and 
' V2? —3 
(2 —1) o/22—7 
2. Differentiate with respect to x 
lo 9/84 @x+aer+23 
6 at+ar+a2 ~ 
3. Find the differential coefficient with regard to 2? of 
n/ x2 —2 


x 
4. Differentiate 


a\%, c . z 
(1) a (=) log A with respect'to log= ; 


(2) log — , with respect to 2, 
(3) x(os2? with respect to x. 
- dy 
5. If (1-—«)¥+(1—y)*=0, find ria 
6. Obtain the differential coefficients of 


(1) log = = (2) e**, with respect to x. 
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CHAPTER XII 


1. Find (1) the nth differential coefficient of 1) with respect to 


1 
@+1) @ 
x, and (2) the second differential coefficient of log, (1+) with respect to log, x. 

2. Find the zth differential coefficient of e* (x — 2). 


3. Having given that 


Fd U +y=0, 
n+l 
prove that a pa =+(an+l)#.9 Tait + (nt41) F9 J_0, 


4. Obtain the second differential coefficient of log (a +bx+ cx?) with respect 
to Vz. 


5. Find the nth differential coefficient with respect to x of 


2272-52743" 


CHAPTER XIII 


x zw. xz? gf 


1. Prove that = ae 3+ i537 790°" 
and show that no odd power of z beyond the first can occur. 
2. Prove that Wile-Aie 23 all eo a—gatt oe eee 
Ji+22 


1 
3. Prove that the first three terms of the expansion of (1+.)* in powers of 


xz are ¢— = + ther. 


4, Expand log (+a? +2) in ascending powers of x. 


5. Expand oS in ascending powers of z as far as the term in- 
volving 2°, 
CHAPTER XIV 
1. Evaluate ut | ers" if 
“> 


2. Find the maximum value of x (x —1) (« —2) between the limits 0 and 1. 


3. Find the values of x at the points where the graph of the function 
(1+?) e!-* has its greatest slope. 
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4. A man in a boat at sea, 5 miles distant from the nearest point of a 
straight shore, wishes to reach a place 12 miles distant along the shore, 
measuring from this nearest point. At what point should he land to reach this 
place in the minimum time, if he can row at 3 miles an hour and walk at 
4 miles an hour? 


5. Given /(0)=1876, f(1)=777, f(3)=19, and / (6) = — 218, interpolate the 
values of f(2), f(4), and /(5) and find the values of x for which f(x) is a 
maximum or minimum. 


6. Find the minimum and maximum values of 
xt — 473 — 87? + 48x — 48. 


7. A window is in shape a rectangle with a semicircle covering the top. If 
the perimeter of the window be a fixed length J, find what is its maximum area. 


z —~a+vagm_— qm 


/ or — q@ 


1 
(2) Lt (1+25)*, 
z—>0 ? 


8. Find (1) Lt 
tra 


CHAPTER XV 


3 
1. Show that sd = A'/(2— #8) approximately. 


By considering the function f(z)=a+bxr+c and using the above relation 
prove that log,c=ct—c™* approximately, where c is a small quantity. 


AA & 
2. Show that A0™— 5 a 0” ...=0, when m> 1. 


3. Prove that, if f(#) be a function the fourth differences of which are 
constant, 
df (x) 


LS! =e LF e-2)- Bf (@-1) 4 Bf (w+) -f(e+2)] 
and hence find an approximate value for es [log f(x)] where 


f(a—2)=42°699, f(w—1)=40°365, f(x)=37°977, 
f(e+1)=35:548, f(w+2)=33-075. 


CHAPTER XVIII 
1. Find the integrals 


a x og (log 2) 
(1) ls dx. (2) me ad (3) fear. 


2. Evaluate | wiatdx. 


H. T. B.I. 10 
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3. Evaluate 
PAL log xdx 
G) Ibid: © = Ce en ers 
' 274? —1712+ 256 
4. Find the value of Ree ada. 
5. Evaluate 
sie dx ; | log (1-2?) dx ; | a” (log x8 dx. 


CHAPTER XIX 


1. If f(z) =a+bx+cx*, find expressions for f(3), se (0) and I, Ft (#) dz in 
terms of f (0), f(1) and f (2). 


2, Evaluat eet 
ls valuate f (Qx+1) RH 
; 1 272492 -—2 

3. Find the value of | 1 FG-e@al) dz. : 2. 
4, Find the average value of x (x—1) (x —2) between the limits 0 and 1. 

1 
5. Find the value of \ Bae: 

0 

CHAPTER XX 


1. From the table 
f(0)=217, f(1)=140, £(2)=23, 7(3)=-6, 


find an approximate value for , f(«) dx, and explain why the result differs 
from f(0)+f (1) +f(2)+f (3). 
2. Prove that approximately 
[7 F@ aed {f(-1)+22f0)4F 0), 
and find thereby an approximate value for 
| . 10,000 x 4—* x 23" dx. 


3. Discuss the error in assuming that 
[f@ 4-7-2144 0) 
if (1) f(x) =a+ba+cx%, (2) f()=100x 4%, 
4. Given f(—2), f(0), £(2), find an approximate value for 
fo soves 


¢ 
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5. Prove that 


i f (#) da= py [5f (1) +8f(0)—f(—1)] approximately. 


If the speed of a train on a non-stop run is as shown in the table below, find 
the approximate mileage travelled between 12.0 and 12.30, using the above 
formula. 


Speed in 
miles per hour 


24°2 
35°0 
41°3 
42°8 
39°2 


6. If f(1)=4157, f(2)=4527, f(4)=5435, find approximately the value of 
4 
| FE 
1 


CHAPTER XXI 


1. » persons are sitting at a round table, and from them three are selected 
at random; show that the chance that no two of those selected are sitting next 
(n—4)(n—5) 

(a—1)(n—2)" 


2. A heap of playing cards contains 6 hearts, 5 spades and 4 clubs. A card is 
chosen at random 9 times in succession and is not replaced. Find the chance 
(1) that there are no hearts left, and (2) only hearts are left. 


3. A man has two sovereigns and four shillings in his pocket; from these 
he selects at random 3 coins, and again chooses from these 3 coins at random, 
the first coin for A, the second for B, and the third for C. 

Find the values of the expectations of A, B, and C@. 


one another is 


4. Find the probability that out of 5 persons aged 45 exactly 3 will die in 
a& year, 
Probability that a person aged 45 will die in a year="01224, 


5. A die with six faces is thrown three times and the sum of the throws is 
twelve: find the chance (1) that the first throw was a four: (2) that four was 
thrown each time. 


6. Two men, A and B, each draw a card from a well-shuffled pack of 
playing cards, find that they are of the same value, and replace the cards; 
they do this four times in succession. Find the chance that this would happen 


and show that it is approximately sx006" 
7 


\/ 


BeatriceGloria_personal library 
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7. The probability that a man aged 50 will survive one year is ‘98428. 
Show that the probability that, out of 5 men aged 50, 3 at least will i within 
a year is 0000385. ‘ 


8. The 26 letters of the alphabet are placed in a bag. A and B alternately 
draw a letter from the bag, the letters drawn not being replaced. The winner 
is the one who draws most vowels. A starts and draws a vowel with his first 
draw. What is his chance of winning? 


9. Ifa number of five figures containing any five of the ten digits once only 
is written down at random, what is the probability that it is divisible by 9? 


10. Given the following table find the probability that one at least of three 
persons aged respectively 20, 30, 40 will die between the 10th and 20th year 
from now: 


Probability 
of surviving 
10 years 


Probability 
of surviving 
20 years 


Age 


11. Two persons, A and B, play for a stake, each throwing alternately two 
dice, A commencing. A wins if he throws six, B if he throws seven, the game 
ceasing as soon as either event happens. What ratio will A’s chance of winning 
bear to B’s? 


12. Thesum of two positive integers (excluding zero) is 100; find the chance 
that their product exceeds 1200. - 


13. The following table shows the probability that a woman of the age 
yn will marry in a year: 


Age Probability of marriage 
20 0665 

25 1033 

30 0649 

40 0183 


Find the probability that, out of 4 women aged 20, 25, 30, 40 respectively, only 
one marries within a year. 


14. A bag contains 8 counters, numbered 1 to 8. Four are drawn at 
random. Find the chances that 
(1) The sum of the numbers on the four counters amounts to at least 17. 
(2) The counters numbered 2 and 3 are among the four. 
(3) The four counters contain at least two of the three counters numbered 
3, 5 and 7. 


15. A penny is tossed six times. Find the chance that neither heads nor 
tails have occurred three times in succession. 
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16. An urn contains counters marked with the digits 6, 7,8 and 9; and the 
number of times each digit occurs is equal to the value of the digit. If counters 
are drawn one at a time, each counter being replaced when drawn, what is the 
probability 

(1) that the digit 6 is drawn before the digit 9; 
(2) that the sum of the first three digits drawn is exactly 20? 


17. A and B playa set of games, to be won by the player who first wins four 
games, with the condition that if they each win three they are to play the best 
of three to decide the set. .A’s chance’of winning a single game is to B’s as 
2to1. Find their respective chances of winning the set. 


18. The probability of any one of 10 men each aged 30 surviving a year is 
99229. 

Show that the probability that exactly 5 men out, of the 10 survive a year 
is 6°6 x 10-9, 

Find also the probability that of the 10 men one particular man will die 
first and another particular man last. 


19. A point is taken at random within the area bounded by the curve 
y= log x, the x axis, and the ordinates at the points x=1, and x=4. 

Find the probability that the distance of the point from the y axis is less 
than 2. 


20. In a game of whist the dealer found that on turning up the last card he 
had the Ace, King, Queen, Knave, Ten, and 3 other trumps in his hand. Find 
the chance that this would occur. 


21. A and B cut a pack of cards, the player who wins the cut six times to 
be the winner. A, having won four times to B’s once, cuts a five. Find the 
chance that A will be the winner. 


22. Ina line AB of length 3a, a point P is taken at random and then in 
AP a point Q is taken at random. What is the probability that PQ. exceeds a? 


ow bp 
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CHAPTER III 


. ©=12 or —16. 


. ath term=n?+n-—3+3". Sum of n terms= 


cane acne Dale 
6 g™ 
20 
51224, a = 99) Sl oe 


2 °°" "(82 4+-1)(82+4)...(82+ 70)’ 


b et. 2!. 6. 2-2 (m+1)(m+4). 


CHAPTER IV 


» f (8)=75857°2 ; f(4)=77623'8; f(5)=79432:0. 2. 02017. 
. 69215. 4 12,14; 26 42; 138,266. 5 97357. 6. 97453. 
. Af(0)= —1°753875, A?f(0)= —°07725 and A*f(0)=-00125 from which 


the series may be completed. 


. °1954. 
. Af(1)=1°1006, A2/(1)=-1083 and A?f(1)=-0137 from which the values 


of f (x) required may be found. 


. (5) =3214; f(6)= 66%. 


24 — 9234. 372-—do+5 or 24. 2%*—872-182—19. 
log 45=1'65321, log46=1°66276, log 47=1°67210, 
log 48 = 1°68125, log 49=1-69020. 


. 5884346. 


CHAPTER V 
2290°0. 2. 95384. 3. -3165. 4. 01625. 


. The series is (25) 568, £(26)585, (27) 602, (28) 621, 


f(29) 644,  f(30)671, (31) 707, (32) 747, 
f(33)'793, (34) 843, (35) 897. 


CHAPTER VI 


. 2751. 2. 24. 3. 169. 4, 37-2. 


CHAPTER VII 


ae _ 
l-r* n(n+1)(n4+2)... (n+r—2)’ 144° 


CHAPTER VIII 


. *6895. 2. 460. 3. ‘500. 
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CHAPTER IX 
1. (i) 10°389; (ii) 10-475. 2. 3-297, 


CHAPTER XI 
— 44-62? +42 -—21 


1] 
igh ee 
za (1 —log 2) (ei at_ayh ae 


x? (3a* + 2ax + x) 1 


7 3 (a+2) (a? +27) (a? +axr +27) ° * wf x2—2- 
4, (1) a (=) x loga log = + ax log = log= +1 +1 
a a a a , 
1 1 3a (logx)? (log x)? 
©) EBT Oa 


5 —¥(L—x)"1—(1-y)* log (1-y) 


1 
az (1—y)*1—(1—2) log (1-2) * 6. ye a3 ex x (log 2 +1). 


CHAPTER XII 


(- Liha! n! 1 2 (n+1) x logge 
Lo) SE Leas ap ees} © Geer 
2. 6%, qn-2 ees 2)? + Qna (# —2)-+n(n—1)]. 
g Lab-+x (6ca— b?) — bex? — — 2c?a*) 


“ + ba + cx)? 


5.2” 2 
5 (nll genaen ~ warm | 


CHAPTER XIII 


a 325 fo) .38 oni(@r—1)g**! 
a4 it oe A(— Ir orale 


oe 5a? 20975 


TN Tg © at eG 
f CHAPTER XIV 
iat. 2: ge 3, a 3 . 4 4. miles from the nearest point 
axle = Ne VT 
5. f(2)=218, f(4)=0, f(5)= -19 _ 
Maximum value when z= ok vi0 : 
Minimum value when ot 


6. Minimum values— 3, —128. Maximum value 0. 


2 Mm—N « 
i Or +4)" 8. (1) 4 /= —, ge; @) I. 
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CHAPTER XV 
3. —-06353 
CHAPTER XVIII 
1. (1) YE" (out + 622+80416)4+05 (2) arg log (pas) +63 


(3) fees (log x) —1]+c. 

| 43 1222 247 24 

” log,a ~ dog.ay? * (log, a) as (log, a)! “ (loge = |+e 
3. (1) $[(e+a) V2? +2ax—a? log {(a+a) +2? +202} +c]; 


(2) } (a? — x2)? — a?/a?— az? +c; (3) TBF + log (1-2) +0. 


—3 
. Tena ttle ssi te 


—7\8 
5. log = 1) totes (1 +2) log (1+2) —(1—2) log (1—x)—2r+e; 
at 7 [{(n+ 1) log x}—3 {(n+1) log x} +6 (n+ 1) log x — 6]+c. 
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1. §[8f0)+6f(1)-F(2)]; $[-3f00) +47 (1) -F(2)]; 
rz [5f(0)+8f (1)—f(2)]- 


2 ss. 3, log,6. 4. }. 5. ie 
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CHAPTER XXI 


2. (1) rye; (2) sdys- 3. 7s. 4d., 78. 4d., 7s. 4d. 4. 0000179. 
5. (1) 8s (2) ay C iatex: 8 hit. 
9. $. 10. 31466. 11. 30:31. 12. $%. 
13. -2109. 14. (1) g83 (2) #5 (3) 4- 15. 42. 
16. (1) 2; (2) rds. 17. A’s chance 4§822; B’s chance 98,. 
, 8 log2—3 101,232 
18. 01111. 19. 64 ]og2— 165 20. 309,966,275 ° 
21. 744. 22. 2-4 log3. 
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